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| iscussion of Current-Sheet Approximations in Reference 
to High-Frequency Magnetic Measurements 


Bohdan Kostyshyn and Peter H. Haas 


The proble m of violations of thin, uniform current-sheet assumptiorfS in toroidal coils 
of rectangular cross section is investigated from the viewpoint of high-frequency magnetic 


measurements 


The errors resulting by neglecting leakage flux and by assumptions as to the 


location of effective current sheets are investigated experimentally as a function of permea- 


bilitv, core dimensions, 


1. Definition of Symbols 


u= Magnetic permeability of the material 

this paper, the initial permeability 

Measured inductance of sample coil containing 
core 

Calculated value of toroidal coil inductance on 
basis of uniform thin current-sheet assumptions 

Permeability of material as inferred from induct 
ance ratio of measured to calculated inductance 

Outside diameter, inside diameter, and height of 
toroidal core, respectively 

Value of wu, reached as the number of turns becomes 
very large 

Absolute permeability error. 

Components of Aw, see text. 

Distance of the effective current sheet from the face 
of the toroid 

Correction term to the 


as used in 


measured inductance for 
pu, to equal ue 

Correction term to the 
Me to equal vy 

Mean spacing between adjacent wires on toroidal 
coil 

Correction term of the type Als as calculated from 
theoretical considerations by Rosa. 


calculated inductance for 


2. Introduction 


It has commonly been the practice in the past [1] 
to determine the initial permeability of a magnetic 
material at radio frequencies by winding a toroidal 
sample with N turns of wire, measuring its induct- 
ance, L,, at low magnetic field strengths, calculating 
the geometrical inductance, Ls, of an air toroid of 
the same dimensions, and taking the ratio of the 
two, L,,/Ls, to be the permeability. Another method 
employed [2] has been to substitute for Ls the meas- 
ured value of inductance (L,,), of an equivalent air 
toroid wound on a polystyrene core of identical 
dimensions, with the same wire and number of turns 
as the magnetic core to be measured. 

The assumption in applying the first method is 
that the given ratio does not depend on the frequency 
(if the permeability is truly constant), the number 
of turns wound on the sample, wire size, and the 
dimensions of the sample core. The assumption in 
applying the second method is that any corrections 
that appear for the inductances for the coils con- 
taining the magnetic core and the equivalent air 
core, respectively, are equal and appear in such a 
manner as to cancel when the ratio is taken. It is 


gures in brackets indicate the literature references at the end of this paper 


wire spacing, and wire size 
quacy of at least two attempts at calculations of these errors 
as to the limitations of measurements on toroidal coils by high-frequency 


Conclusions are reached as to the ade- 
k inally s conclusions are draw n 
considerations 


the purpose of this paper to test the validity of these 
assumptions 


3. Limitations on Test Parameters 


For the test the various parameters were limited 
such that the following conditions were adhered to: 

1. The magnetic field strength was sufficiently low 
that the permeability «4 measured would be equal to 
the initial permeability; that is, that value of the 
permeability obtained as the magnetic induction, B, 
and the magnetic field intensity, /7, simultaneously 
approach zero. 

2. The frequency was low enough that capacitance 
effects from wire proximities could be demonstrated 
to be negligible, yet it was high enough that the 
internal inductance of the wire could be neglected. 

3. The magnetic materials were chosen so that the 
permeability showed no change as a function of fre- 
quency until a frequency many times the test fre- 
quency was reached 

4. Realistic wire sizes allowing a tight, uniform 
winding were used { toroidal winding machine 
was used (fig. 1) to apply single-layer windings to 
cores of rectangular cross section 


Figure 1. Typical coils used, 
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4. Analysis of Assumptions 


We now define the effective permeability of a 


toroidal sample core as 
(1) 


where yu, is effective permeability; L,, is the measured 
inductance of sample coil, and Ls is the calculated 
inductance of an air toroid, equivalent to sample 
coil, assuming a thin, uniform current sheet, where 
Ls Is equal to 


b 
Ls=2N°*h log, 10~* henry, 
a 


where V is the number of turns of wire on sample, 
h is the height of sample in meters, 6 is the outer 
diameter of sample, and a@ is the inner diameter of 
sample. 

If one now varies the number of turns on a given 
magnetic core, with given wire size and plots the 
effective permeability yu, as a function of N, it is 
observed (see fig. 2) that the assumption that the 
ratio L,,/Lsis not dependent onN, as assumed earlier, 
does not hold, but that for small N a large value of 
uw, results, and as N is increased, yu, approaches 
asymptotically a value yw, slightly higher than the 
true permeability u, as determined by an NBS radio- 
frequency permeameter [3] 

One now recognizes three values for permeability: 

uw equals true initial permeability; that value of 
permeability obtained by using eq (1), provided a 
thin, uniform current sheet exists adjacent to the 
toroidal sample 

Me equals effective permeability : that value of per- 
meability obtained from eq (1) when the sample core 
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Fieaure 2. Variation of u, with number of turns, N. 
Wire: Insulated No. 38 Core: Outside diameter, 1.507 in.; inside diameter, 
1.004 in.; height, 0.375 in 


Schematic cross section of the surface of « 
4 1 


coil wound with round wire of diameter d 


is wound with N turns of actual wire of finit 
mensions 
Me equals asymptotic permeability; that val 
effective permeability reached as the number of | 
N becomes very large (fig. 2 
We further define 


Au=Api+Apo= yep, 
Aui=a— bs, 
Aue Me Ma’ 


As noted before, the assumption made in applvir 
eq (1) and assuming it to vield the true permeabilit 
is that a constant current sheet exists about and in 
mediately adjacent to the toroidal core. This cor 
dition requires that there be no external magnet 
fields beyond the confines of the current sheet 

For an actual coil wound with wire of finite di- 
mensions and individual turns spaced a finite dis- 
tance apart (fig. 3) it is seen that a thin current sheet 
may not exist, but one may visualize a hypothetica 
current strip in each wire. If one assumes the ex- 
istence of a definite current distribution, one might 
calculate a mean or effective current strip, located 
finite distance z from the face of the core. It 
further noted that as the number of turns is increase: 
the mean wire spacing D) is decreased and the curren 
sheet is more nearly approximated. One may pr 


dict on this basis that for a few turns the predominan' 
cause for the error Ay, is the existence of a series 0! 
current strips, rather than a current sheet. This 


then implies the existence of magnetic fields beyon 
the confines of the core. For a large number o 
turns one may ascribe the error Ay, to be caused b 
the location of an assumed current sheet at son 


finite distance from the faces of the core, as deter- 
mined by the current distribution in the wire rather 


than immediately adjacent to the core. 


5. Discussion of Correction Terms 


It has been demonstrated in the literature [!, 4 
that corrections might be made to account for thes 
more realistic conditions. 

In one case it is assumed that the magnetic fw 
exists only within the confines of the coil, including 
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space not occupied by the core The perme- 


is then calculated from 


-AL, 


henry {1 


2h 


log, b/aX 1077 henry 


‘ 


nductance of au space between current sheet 
ore; A, is the residual area of air space; d, is 
ve magnetic diameter; and z is the distance of 
fective current sheet from the face of the toroid 

1, one may alternatively apply a correction 
in eq (1), so that the air inductance would 
vely be increased to compensate for the larger 


ive coil cross’ section. Equation 2 then 


2N*(h+-: log, ~~ 10-7 henry 


must be noted at this point that the above 

ections account only for effects of finite wire 
The fact that leakage fields do exist 
ond the confines of the coil is not considered. A 
oretical treatment from the literature of the 
ts of stray fields will be discussed in a later 
tion, after a consideration of additional experi- 
ntal results. 


ensions 


6. Experimental Results 


We define the correction terms AL,, and Als, such 


Le 
” Lsgt+ALsg 


No geometrical meaning is affixed to either correction 
rm, except that when either is applied, all necessary 
correction will have been taken into account. 
It may be shown from eq (8) and (9) that 


AL,,=Apls 


Au L;* 


AL L,»—AuLs 


AL. 
v3 


AL; (12) 
ich relate the correction terms to measurable 
intities. 

if " " 

if now the number of turns on a group of cores of 
same size but different permeabilities is varied 
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and the effective permeability, u,, is plotted as a 
function of N? (fig. 4), it is observed that within the 
limits of the test, Au is the same for all samples for 
any particular value of N. Thus for a given core 
size and wire size, Ay is only a function of the number 
of turns on the core, i. e., 


Au 


With this result in mind, a detailed study of the 
behavior of Aw is made on a polystyrene core, and 
it is assumed that the conclusions that will be drawn 
will be applicable to all other cores of the same 
physical dimensions wound with the same wire 
The reasons for choosing a polystyrene core for this 
detailed study are 

l. w is very accurately known (u=1)[5] 

2. The percentage error Ay/u is largest for the 
polystyrene core, and may therefore be more pre- 
cisely determined 

Figure 5 shows the results of another series of 
measurements on a polystyrene core, similar to 
those for the magnetic cores, with the effective 
permeability, u,, plotted as a function of N, with 
a given core and wire size. It is noted that, indeed, 
the general shape of the curve is the same as in figures 
2 and 4, and that the value of the asymptotic perme- 
ability, u., is 6 percent above the true permeability. 

If one now assumes a current distribution such 
that the current is concentrated on the periphery of 
the wire, and from this deduces the distance z be- 
tween the current sheet and core face to be half the 
diameter of the wire, one may apply the correction 
terms in eq (7) or (2a), substituting d/2 for x. It 
is observed that both correction terms vield approx- 
imately the same result, i. e., the entire curve is 
shifted down and approaches the value of u,—1.03 
asymptotically. Therefore, this new assumption 
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FIGURE 4 
rameter fora group of magnetic cores of given dimensions with 


Variation of u, as a function of N? with was a pa- 


given wire size, 


Dotted lines indicate values of true permeability Wire: Insulated No. 38 








decreases the error to 3 percent, but does not account 
entirely for the error Ay. This indicates that the 
current distribution in the wire that has now been 
assumed might be incorrect because of the proxim- 
ity of neighboring turns of wire, which certainly 
will influence the glistribution. Furthermore, the 
experimental evidence indicates that the effective 
current sheet is located at a still greater distance 
from the face of the core. It might be noted at 
this point that the problem of winding a wire ad- 
jacent to the core face is a difficult one, and it can 
be shown that even a small air space will lead to a 
considerable error Ax. 

Figure 6 shows a plot of », as a function of N, 
with wire size as a parameter. It is observed that 
the basic curve is shifted upward and to the left for 
increasing wire size. The curves indicate an in- 
crease in the value of the asymptotic permeability 
uw, With increasing wire size, showing that the 
current distribution is altered in such a manner as 
to increase the distance of the effective current sheet 
for larger wire. At the same time, because of a 
greater effective “current strip width” per turn, the 
uniform current sheet is more closely approximated 
by the heavier wire for small values of N. This 
causes less flux to be lost from the test core by leak- 
age, which in turn results in a smaller Ay. Thus 
it may be seen that for small ', the smaller wire 
yields a larger Av, whereas for large N, the smaller 
wire yields the smaller Au 
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Fieure 5. Variation of u, as a function of number of turns N 
for a polystyrene core. 


| 
Wire: Insulated No. 38. Core: Outside diameter, 1.504 in.; inside diameter, | 
1.004 in.; height, 0.373 in 


If one plots u, as a function of 1/N for a 
styrene core (as in fig. 5), it is observed that 
limits a linear relationship exists. If, further 
the core size is varied (see table 1), and y, is p 
as a function of 1/N, with a given wire size ({ 
it is observed that the slope is the same for all 
tested within the linear portion. Considering 
the linear portion of the curve and choosing thy 
permeability as a base point, it is concluded 
within the limits of the test, Ay is a linear fw 
of 1/N, 1. e., 

ae , 
Au N + C\. 


Figure 8 shows a plot of u, as a function of 1/A 


wire size as a parameter; it is observed that 


slope Kk, is dependent on wire size. 
K f(d). 
Dimensions of the cores used in jigu é 


TABLE 1 
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Ficure 6. Variation of pu, as a function of number of turns N, 


with wire size as a parameter. 


Core: Outside diameter, 1.504 in.; inside diameter, 1.004 in.; height, 0.371 in 











The term C, is numerically small by comparison to 
A, and apparently dependent on both core and wire 
\ mean spacing between wire enters may also be 
fined by 
T b a 


D aN 16 


\s D) is linearly related to NV, it may be concluded 
it within the range stated, the error is a linear 
ction of the mean spacing between wire centers 
is eq (14) may be written 


Au=K,D+C,, 


9) A 
K, 2K . (18) 


r(b ra 


data of figure 7 is replotted on figure 9, which 
ws a plot of u, as a function of PD, with core size 


ws 


Vi a | 
Pd no~ Ya 
YO " 








as the parameter. It is interesting to note that Ag 
is largest for the smallest test sample, with a given 


wire spacing, D 


7. Effects on High-Frequency Magnetic 
Measurements 


Sufficient experimental evidence has now been 
presented to allow us to draw some conclusions as 
to the effect of flux leakage on permeability measure- 
ments at high frequencies. The quantity of major 
interest in an analysis of the accuracy of a magnetic 
measurement is the fractional permeability error, 
E, given by 
Au 
M 
From eq (14) we have then 
ky , ¢ 
= 20) 
uN p 
The values of A, and C, obtained for a number of 
different core and wire sizes were all such that for 
the number of turns NV usually employed at frequen- 
cies in the megacycle region, the error can be approxi- 
mated by 


Au kK 


m4 ~E 21) 
u pN 
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Under these conditions, one then obtains the condi- 
tion for constant error 


uN constant 


for a coil with given wire size 

Let it now be assumed that a toroidal core of per- 
meability 4, has been wound with N;, turns, yielding 
a measured inductance L, Let it furthermore be 
assumed that an error /, (this error caused by viola- 
tion of current-sheet assumption only) results from 
a calculation of permeability by using eq (1), and 
that this represents a maximum error that should 
not be exceeded in the measurement of other cores 
of permeability u, of identical dimensions 

From eq (8) we have 


ui bs. +ALL 


From eq (10) and 
! 
(24) 


where / Au,/m, and Ls=N?°L Now a core of 
permeability uw, is to be measured with an error, F), 
no greater than /,. If we consider the error to be 
held constant, i. e., A= &,, then from eq (22), we 
have 
u . 
N 
Me 
We now investigate the value of L,,.. i. e., the 


measured inductance of coil 2, to satisfy the constant- 
error condition. As in eq (24), we have 


(1+E, 


solving for 


Dividing eq 26 by eq 24 


Len 


—L 


ba 


This is the value of measured inductance necessary 
in sample 2 for the error Ay:/u, to equal the error 
Au, Mi. 

If permeability My of core 2 is assumed to be less 
than that of core 1, it is noted that for the errors 
from violation of uniform, thin current-sheet as- 
sumptions to be the same, coil 2 must possess an 
inductance larger than that of coil 1 by the ratio of 
the two permeabilities. 

Up to this point only errors in the determination 
of high-frequency permeability caused by violating 
current-sheet assumptions have been considered. 
Errors, however, are also introduced by the presence 
of capacitance between the windings of the coil. In 
particular, if the inductance, L, 


of a coil, and its | 


distributed capacitance, C, may be regard 
lumped inductance and capacitance in parall 
effective inductance of the coil is given by 


L 
L. wbC 


l 
2.yLC 


where L, is effective inductance; LZ is actual in 
ance; fo is the self-resonant frequency of th 

In order for the measurement to remain reaso1 
free from errors caused by distributed capacita 
as shown in eq (29), the measurement frequen 
must be much lower than the natural resonant 
frequency, fo. 

Considering the case discussed above, wherein th 
error Au/u for cores of permeability yu; and yp is to by 
constant, uw, being less than y;, it is noted that coil 2 
must have a higher inductance than coil 1. How- 
ever, in order to have this larger inductance, mor 
turns had to be placed on core 2 than on cor 
resulting in an increase of capacitance for coil 2 
Consequently, coil 2 will have a lower natural reson- 
ance frequency fy than coil 1 because of both highe: 
inductance and capacitance. It thus appears that a 
smaller measurement error will result for a higher 
permeability at a given frequency. 

It has been shown by numerous authors [2] that 
the effects of distributed capacitance may be cor 
rected for by measuring the self-inductance at several 
frequencies. It must be pointed out, however, that 
this procedure assumes a constant self-inductane 
and thus a constant permeability with frequency 
If this, however, be known a priori, a meausrement 
at a frequency where errors from residual capacitance: 
may be neglected would suffice. 


8. Experimental Evaluation of Theoretical 
Corrections 


As noted previously, attempts have been made t 
correct the approximate eq (1) for the assumptions 
made 

In the formulation already described, it is inferre: 
that there are no leakage fields existing beyond th: 
confines of the coil. If this were so, then the cor- 
rection AL, in eq (6) would equal the correctio1 
AL,,. From eq (7), we see that 


AL,.=CN’. 
It has been shown in eq (14) that 


ky 


Au y tO. 


But we also have from eq (10) and from the iden- 
tity Ls=N’L,, 


AL, =ApN’*L,. 
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this equation and eq 14), we have 


AL, =K;N+ K.N’. 32) 
all but large values of N, AL,, is dominated by 
rm K,N, and only for very large N is AL, 
nearly a function of N*. Therefore, it may 

meluded that eq (30) will apply only at large 
here the assumption that a uniform current 
exists is more nearly valid. However, for high- 
ency measurements it is important, as stated, to 
the distributed capacitance down to a minumum 
which is seldom feasible if a constant, uniform 
nt sheet is simultaneously required 
other correction based on a theoretical approach 

s into account the possibility of leakage fields 

ting, when the core is wound with round wire of 

Under such conditions the self- 

ictance is affected by those lines of flux within 

wire and by those linked with each separate 

of wire in addition to those lines of flux within 

test The correction worked out by Rosa 

rs the calculated value of inductance L,. The 

rection term depends on the pitch and size of 

and may be positive or negative. It is deter- 

ned by considering the difference in self-inductance 

mutual inductance between turns of a coil of 

round wire of finite dimensions from those of a coil 
of thin tape. 

Rosa’s equation for the self-inductance of an air 

toroid is [4] 


dimensions. 


core 


—T 
A log. (1 745 D , 
where / is the length of one turn of wire equals 
2h+b d is the diameter of the wire used; D is 
the mean wire spacing (see eq 16); B is the correction 
from Grover’s tables [6] for mutual inductance be- 
tween turns We define Rosa’s correction term 


AL 2NI(A 

lt is noted at this point that if AL, included all of the 

necessary corrections, it should equal the correction 

term AL, for an air toroid, (see eq 9 We thus 

onsider a polystyrene core and from the observed 
ta, plot AL, as given by eq (11 


Appl? 


, (11 
—~Auls 


ALs _ 


me uses L, as defined by eq (2) to determine 

a plot of AL, (see top curve of fig. 10) as a fune- 
of N does not compare favorably with a plot of 
Rosa correction term AL,. Furthermore, if we 
stitute a corrected value of Ls, assuming the 
ent sheet to be a distance d/2 from the sample 
face, it is noted that the agreement is but slightly 
roved. Recalling our earlier prediction as to 
iration of permeability errors into two parts, and 
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Fiaure 10. Compa 


on the basis that AL, takes care of nonuniformity of 
current sheets only, we calculate ALs by using 


Aw as shown: 


(lla) 


We now plot this equation as a function of V (see 
fig. 11) and notice that a favorable comparison with 
Rosa’s correction term AL,» is obtained. 


As a result of eq (10) and (12) and the experimental 
fact that Au is not a function of u, we find that 


Auls 


vi 


(36) 


ALs 


Thus the values of AL, for two cores of permeabilities, 
My and Me, are related by the equation 


ALs.=*" 
- Me 


ALs. 

While Rosa’s correction term is strictly valid only 
for an air toroid, it was pointed out earlier that in the 
calculation for an air toroid it takes the place of our 


correction term Als. In comparing its possible 
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validity for a coil with arbitrary permeability «, we 
thus let it take the place of Als, in eq (37), with ug 
now equal to unity. The two curves, that is, 
AL, for a magnetic core and AL, divided by that 
core’s permeability, are shown plotted versus .V as 
the bottom curves on figure 11. The fairly good 
agreement obtained indicates a possibility of esti- 
mating a value of the permeability for a magnetic 
core by using Rosa’s correction term in the manner 
indicated 

We have now shown the effects of two possible 
calculations of correction terms. Referring once 
more to eq (4 


Au= Ap, + Ap, 


it may be stated that the error Ay, is indeed associated 
with violation of the assumption that the current 
sheet is adjacent to the toroidal sample. As shown, 
this may be accounted for by eq (7) or (2a), provided 
one knows the distance z from the face of the test 
sample to the location of the effective current sheet. 
The value of Ay, may be associated with the error 
caused by the presence of leakage fields beyond the 
confines of the coil. This error may be estimated, at 
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least for an air coil and within limits for any ec 
applving the Rosa correction term 


9. Permeability Measurement by Air- 
Substitution 


Another method mentioned earlier for meas 
high-frequency permeability utilizes the mea 
inductance of a coil wound on a nonmagnetic 
This coil is to have the same physical dimen 
as the magnetic core whose permeability is | 
and wound with the same numly 
wire For the permeability, one 


determined 
turns and 


takes 


L. 
L. 


However, by referring to eq (8) and (12), we 


Ln )m=u( Ls+AL 


L)(L AL. 


and eq 38) becomes 


u(Ls ALs_) 
Let+tlhe 


Only when ALs =ALs, will this ratio reduce to 4 
On the basis of the experimental evidence and 


10) and (12), one obtains 


AL s. +s 


AL, =AuLs 


becomes. substituting 


and (43 


Equation 4] upon 


results of eq (42 


\ Au u 
we 1+Ap 


As AuxAp/u, unless p=1, pez, will always lead 
values on the low side, the absolute error becoming 
larger for higher permeability. A practical difh- 
culty that may be encountered is that it is extremel) 
difficult to wind coils exactly alike, and thus further 
errors might be introduced. 


10. Experimental Details 


The magnetic cores used for this investigation 
were of powdered iron. Field strengths abov 
oersted are usually required to effectively change 
the permeability of such cores; the field strengths 
emploved were below 0.01 oersted in all cases. Car 
was taken to assure that changes in the field in- 
tensity by factors as large as 2 to 1 would not chang 
the coil inductance by even as much as | part In 
10,000 





frequency used was 50 ke. The coil with 
ehest reactance used in this investigation 
easured in the range from 1 to 50 ke The 
inductance showed a change of less than | 


LO.000, W hich could be ascribed to capacitance 


Ve 
bridges used for the measurements were two 
Maxwell-type inductance bridges de- 


juality 
Connections 


at Bell Telephone Laboratories 
en the unknown terminals of the bridges and 
ils were made in such a manner that proximity 
were held to a 
lead inductance 


neighboring objects 
that the 


s from 
ble minimum 
be accounted for 

true initial permeabilities of the magnetic 
determined by a_ specially calibrated 
permeameter used in conjunction 


and 


were 


Trequeney 
alibrated capacitors 

order to apply the windings as uniform and 
as possible without any crossovers whatever 


no visible injury to the insulation, a toroidal 
winder capable of applving this type of single- 
winding was emploved Figure 1 shows an 


iple of such windings 


WASHINGTON, January 7, 1954 


11. Summary 


The errors encountered in high-irequency magnetic 
measurements have been discussed mostly from the 
viewpoint of violation of uniform, thin current-sheet 
assumptions. The absolute error in permeability 
has been shown to be independent of the permeabil- 
ity of the material tested, and its dependence on the 
spacing between turns has been formulated. The 
adequacy of certain correction terms suggested by 
other authors have been investigated. It has been 
shown that an approximate correction may be made 
terms, pro- 
into space 


combination of several such 
the leakage flux 


utilizing a 
vided the 
estimated 


extension of 
may be 
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Thermodynamic Functions for Carbon Dioxide in the 
Ideal Gas State 


Harold W. Woolley 


The thermodynamic functions C°/R, (H° 
dioxide in the 
data. 


1. Introduction 


hermodynamic functions for carbon dioxide were 
puted by Kassel 1], whose values were corrected 
ore recent values of the physical constants and 
nded from 3,500° to 5,000° K by Wagman, et al 
Values by H. L. Johnston have been given by 
Kulik, and Gay [3]. A subsequent cal- 
ition by Gratch [4] differed appreciably, showing « 
newhat different behavior at elevated temperature. 
\Vith improved spectroscopic constants, it has become 
sible to carry out a recalculation of greater accur- 
Values obtained for the specific heat in the 
iin are intermediate between the previous cal- 
ations 


nwood, 


2. Spectroscopic Constants 


The spectroscopic constants used for C?OY with 
vibrational energy formuls are 3 1342.86 cm 
2.20 em", t= —0.75 cm 
| given by Taylor, 
549.15 em 


667.30 em. 2 
}.76 em, and 4 1.03 em 
Benedict, and Strong [5], with »,-=2: 
-18.54 em!', and 2,;=—12.48 cm, as indicated 
a -private communication from Benedict, and 
-12.46 em™', chosen to fit the observed 11496.5 
at v,=5. The rotational constant By was 
ken as 0.39038 cm received in a private 
nication from Earle K. Plyler, of the Bureau, 
liminary value. For the stretching constant 
),, a theoretical value of 0.132 107° cm~ was used 
the ground that it was not ruled out by a pre- 
ninary empirical value of (0.165 em~'+0.02 em 
0 of Ply ler For rotational-vibrational inter- 
tion, the values a, =0.00056 em™ and a.2= —0.00062 
were given by Herzberg [6], and a,;=0.0031 
was obtained from Plyler’s newer data. The 
rmi resonance constant W, was taken as 51.01 
correcting for a factor publication 
lay lor, Benedict, and Strong. 


com- 


as a 


in the 


3. Method of Calculation 


The calculation was performed by computing the 
rmonic-oscillator rigid-rotator thermodynamic 
netions and adding corrections for deviations 
m this simple model. Corrections were made 

anharmonicity, rotational-vibrational interac- 


s in brackets indicate literature references at the end of this paper 


Ey) / RT, f Ey)/R 


ideal gas state are given from 50° to 5,000° K, | 


T, and S°/R for carbon 
vmased on the latest spectroscopic 


tion, azimuthal quantum effects, rotational stretch- 
ing, and Fermi resonance. The first three of these 
corrections are an extension of a method indicated 
by Maver and Mayer [7] for diatomic molecules and 
used for polyatomic molecules by Stockmayer, 
Kavanagh, and Mickley [8]. The general formulas 
as now extended are given in another paper [9]. 
The constants as given in section 2 are in form for 
this with the exception of », for which the 
derived formulas require a.=»,—2; In the present 
application to carbon dioxide the anharmonicity 
formulas were carried through the second-power 
terms (including cross products) but not through 
the third-power terms. The rotational-vibrational 
interaction was carried through the first power, 
second-power terms being omitted. The azimuthal 
B+-gy or B+x, terms were carried through the 
second power, third-power terms being omitted 
No correction made for the small splitting 
associated with azimuthal /-type doubling, on the 
that its effect is entirely negligible. C 
product terms between azimuthal and rotational- 
vibrational interaction, between azimuthal and 
anharmonic, and between rotational-vibrational in- 
teraction and anharmonic constants were found to 
have small effect below 1,500° K and were omitted 
throughout. It may be supposed that a complete 
and exact representation of the energy levels would 
include corresponding interaction terms, so that the 
lack of such empirical constants may help justify 
the omission of related cross-product terms. In 
addition, the considerable uncertainty in many of 
the available constants provides a greater uncer- 
tainty in the resulting table than the effects of the 
omitted cross-product terms. The procedure used 
in calculating the effect of Fermi resonance and the 
results of that calculation are given in a separate 
publication {10}. No corrections were introduced for 
either vibrational or rotational cutoff 

The correction quantities contributing to the 
thermodynamic functions were calculated for C’O} 
at 300°, 600 1,000°, 1,.200°, 1,500°, 2,000°, 3,000 
4,000°, and 5,000° K, and for the specific heat at 
several lower temperatures. The harmonic 
lator contributions were computed at closer inter- 
vals and adjusted to frequencies representing a 
weighted average according to the relative abun- 
dance of the isotopes. The molecular weight and 
rotational constant were similarly taken for the 
natural isotopic mixture. The values of C°/R 
interpolated to the intervals given in the final table 


use, 


was 


ross- 


basis 


oscil- 
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allowed tabular integration to give the other fune- 
tions, duplicating the directly calculated values to 
about one in the fourth decimal place up to 1,000 
K and satisfactorily at higher temperatures S°/R is 
given as the sum of (// I: RT and I E> 
RT, as tabulated. The translational contribution to 
the latter is (5/2) In 7 3/2) In M+In (22°?) 
(A? N2°P.,), evaluated as (5/2) In 7'+2.0121, 
on \/=44.01 for molecular weight and values for 
the physical constants. The value used for he/k is 
1.43847 em deg 


based 


4. Thermodynamic Functions 


thermodynamic functions calculated for 
dioxide include heat capacity, enthalpy, 
and entropy for the standard ideal gs 
listed in the dimension- 

RT F"—E,)/RT, 


corrections added to the 


The 
carbon 
free energy 

The properties are 
forms (€)/R, H Ek 
in table | The 


state 
less 


and S°’R 


RELI 

0645 
0197 Zh. 457 
O44 23. 1653 
08670 23. Mao 
0905 23. 3500 
1139 23. 4533 


1371 23. 5460 
1602 23, 6372 
1831 23. 7280 
ys) 23. 81452 
2285 23. #022 


2h 23. GR7R 
7 24. 0722 
24. 1554 
24. 2374 
24. 3183 


harmonic-oscillator rigid-rotator functions bas 
vy, and »;, to obtain C°/R, (H°—E°)/R7 

F°— B)/RT for C’O}”, are given in table 

The amounts by which the functions fo 
natural mixture (o nitting entropy of isotope 1 
and of nuclear spin exceed those for C “O,"" are 
in table 3. 

Although it was desired to produce an ace 
table of thermodynamic functions, the new \ 
considerable uncertainty. Using 
trarily assigned uncertainties for some of th. 
stants, rough estimates of the magnitude o 
uncertainties for the functions C°/R, (H°—E 
and I k.)/RT are listed in table 4 for 
1.000 3,000°, and 5,000° K. Uncertaintir 
0.1 em for » and 2» have been taken arbitr 
and may be too small. <A 2-percent uncertaint 


have some 


the Fermi effect as computed is consistent with t) 
order of magnitude of 
values of W 


differen: 
Her: t 


the 1.3—em 


different For a and a 





Values of coefficients of change of thermodynami 
properties with respect to small changes of spectro- 
scopic constants are given in tables 5, 6, 
for the anharmonicity, rotation-vibration interaction, 
and the azimuthal B+z,, constants. If 
revisions of these constants involve 


and , p 


ensuing 
only small 
changes, the corresponding revised thermodynamic 
functions can be adequately approximated with 
small proportionate changes. 
cates an uncertainty of +0.00010 cm™', or about 
sixth of their values. For each temperature 
sum of the separate uncertainties has been 
tained It is to be recognized that the actual 
ertainty is probably considerably greater in view 
the arbitrary values assigned and the neglect of 
of the constants in making these estimates. 
\ comparison of table values with the experi- 
ntal values reported by Masi and Petkof follows 


CeR 
» 


experimen tal 








Tape 6 Coe ficients for change of 


H E; PT with change 


of spectroscopic constants 


Values are 10° ¢d{(11 Ff 


hk 

00 0, O01 

ow 10s 
1,000 .4 
1,20 1.00 
1, 5001, 57 y 
2,000 2.5 13 
4,000 4.5 23 
ioe. 4 2 

000 8.3 41 


Coefficients for change 


change of spectrosco} 


Values giver 


» OS9/0 
70 ool 
ow 


RT )/de 


0. 001 
161 


us 


F SIRT with 


constants 


0. 136.0 

Rl o2 
1. oS 280) 
2.6 Al 


$f us 


5.2 Leo 
a4 38 
1.7 60 
15.1 8.2 


[1] 
[2] 
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Optical Spectroscopic Determination of Hydrogen 
Isotopes in Aqueous Mixtures 


Herbert P. Broida, Harold J. Morowitz,' and Margaret Selgin 


Extension of the optical spectroscopic method of isotope analysis of gaseous mixtures 
of hydrogen and deuterium to mixtures of water and heavy water provides a direct method 


for determining the total water content of materials 


A determination can be made on a 


sample of less than 0.1 milliliter in less than 10 minutes to an accuracy of a few percent 
Using a sample of at least 1 milliliter, an accuracy of the order of 0.1 percent is obtained in 


an hour or less 


This paper describes the system and the investigation of optimum operating 


conditions and effect of variations in operating conditions on precision and accuracy of 


measurement, as well as limitations of the technique 
of distilled 


method can be used to determine the water content of most materials, 


measured directly almost as easily as one 


liquid, or gas phase 


1. Introduction 


Following initial work on an optical spectroscopic 

alysis of gaseous mixtures of hydrogen and deu- 

um [1,2],? a similar technique was tried with 

<tures of water and heavy water. A preliminary 
study [3] indicated that it was possible to carry out 
such analyses to an accuracy of better than 3 percent. 
lhe work discussed in the present paper represents 
in extension of the preliminary study in an effort to 
mprove both the precision and the range and to 
determine present limitations of the analysis. A 
detailed study was made to determine the influence 
of operating variables on precision and accuracy of 
measurement in order to find the optimum operating 
conditions In these studies the following accuracy 
was obtained: (a) 0.1 percent of the intensity ratio, 
in the isotope ratio range 0.01 to 1.0; (b) 1 percent, 
in the range 0.00015 to 0.01 (see table 5) 

In addition to the above the rather practical prob- 
lems of developing the method into a routine analyti- 
cal technique were taken into consideration. For 
this reason efforts were made to achieve a rapid anal- 

sis, and operating conditions were viewed with the 
possibility of constructing a simple flow system and 
optical system for routine use. Several types of 
routine measurements were made concurrently with 
the study of the technique 

lhe technique used is as follows: The vapor from 
an aqueous sample containing a mixture of H,O and 
D.O is pumped through a tube where a high-fre- 
quency electrodeless discharge is maintained. The 
light emitted from the discharge tube contains lines 

the Balmer series of hydrogen and deuterium with 

very small amount of background radiation. The 

emitted light is dispersed by a monochromator, and 

relative intensities of the 8 lines of the Balmer 

ries of hydrogen and deuterium are recorded and 

easured. The ratio of the intensities of the two 

esent address, National Heart Institute, Nationa] Institutes of Health, 
sda, ( 


gures in brackets indicate the literature references at the end of this paper. 
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A sample of blood or serum can be 
water The optical spectroscopic 
whether in the solid, 


lines is a monotonic function of the relative concentra- 
tion of the isotopes, and to a high degree of accuracy 
is equal to the ratio of the isotope concentrations 
For precise work the measurement is a relative, 
rather than an absolute one, in that the relationship 
between intensity ratio and concentration ratio is 
obtained from calibration curves taken on standard 
samples whose compositions are determined by inde- 
pendent measurements 

The statistical terms [4, 5] used in this report are 
defined as follows: ‘‘Precision”’ refers to the variation 
encountered during readings taken in close succession 
under identical conditions. ‘Reproducibility’ re- 
fers to the variation among readings when they are 
spread over a considerable time period, It is not a 
fixed quantity but depends somewhat on the time 
interval and the history of the apparatus in that 
interval. “Accuracy” refers to the variation of the 
measured concentration from the true value and is 
determined from the precision and reproducibility 
of the measurements and from the accuracy of the 
calibrated samples. The relationship between repro- 
ducibility and time determines how frequently it is 
necessary to take calibration points These fae- 
tors are all dependent on the stability of operating 
conditions. 

The above variations are expressed either as stand- 
ard deviation or as the coefficient of variation. These 
terms apply either for single measurements or for 
averages. As the standard deviation of the average 
is dependent on the number of readings taken, this 
value is most useful when compared to other sets of 
the same number of readings. For the purposes of 
this study, all standard deviations of the average 
refer to sets of 10 readings. 

Stability of this system may be viewed in the fol- 
lowing manner. Consider the measured ratio of 
intensities of the D8Balmer line to the H8Balmer 
line at 4860.0 A and 4861.3 A, respectively, and rep- 
resent this ratio by D,/H,. This ratio is a function 
of a large number of variables such as the actual 
concentration ratio in the sample D/H, the tempera- 








ture, the pressure, the flow rate, etc., and may be 


represented in the following fashion: 


D JH,=D/H D/H, X,, X } # 
In 
that 


where 

prac Lice 
f is verv nearly unity. 

, \ completely stable set of m serge r conditions is 
in which D,/H, is a function of D/H only, and 
is unaffected small changes of X,, X, . p a 
From a practical point of view, the 

stability are that changes of D,/H, due 
of the ’s are small compared with changes due to 
variations of D/H. Many the experimental 
curves shown in the following sections are graphs of 
the ratio D,/H, as a function of one of the independ- 
ent variables, with all other variables kept as constant 


Ay are 


operating conditions are 


independent variables 
found such 


one 
by 
conditions of 
to variations 


of 


as possible 

The experiments are 
in all cases aqueous solutions or SUSPCHSIONS In which 
the solvent is a mixture of hydrogen oxide, HO, and 
deut rium oxide D () The number used to describe 
unless specified, the 


D/H. This 


samples employed in these 


concentrations Is, otherwise 
atom of deuterium hydrogen, 
figure depends only upon the number of atoms of 
hvdrogen and deuterium in the sample and is inde- 
pendent of the equilibrium between H,O, HDO, and 
DO 

In taking measurements of ratios 
vapor is pumped from the liquid through the 
charge tube. This procedure raises two questions; 
what is the relationship between D/H in the liquid 
and the vapor and what fractionation (isotope sepa- 
ration) occurs in the flow The first 
tion is purely academic with respect to the analysis 
of water. A given ratio of D to H in the liquid will 
always give rise to a reproducible ratio in the vapor, 
and as all measurements are taken relative to cali- 
bration samples, the measurement taken in all 
cases, referred back to the liquid 

Because of the difference in vapor pressure between 
H,O, HDO, and D,O, the process of evaporating 
liquid from samples involves some degree of frac- 
tional distillation. This has been investigated experi- 
mentally to study the over-all effect and will be 
discussed in a later section. 

The processes occurring in the breakdown of water 
vapor in the radio-frequency field are imperfectly 
understood and will be discussed only briefly. The 
discharge through the water vapor in the tube 
causes excitation, ionization, and dissociation, giving 
rise to the following atoms and molecules from a 
mixture of normal and heavy water: 


ratio to 


the 
dis- 


isotope 


process? ques- 


Is, 


H, D, HD, Hy, D,, OH, OD, H,0,, 


D,O, HDO,, 


D,O,, H,O, HDO, 


as well as other groupings, most of which may exist 
as normal, excited, or ionized molecules. The emis- 
sion spectrum from the discharge contains lines and 
bands from many of these components. The Balmer 
lines of hydrogen and deuterium are prominent in 
this spectrum. 
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The 8 line of the Balmer series was chosen be 
the largest output signal was obtained for this 
This is a result of the product of the actual inte 
and the phototube sensitivity being a maximu 
the region of the 8 Balmer line of the hyd: 
spectrum. The a@ lines at 6561.0 A and 6562 
for D and H, respectively, also can be used; this 
has the advantage of a somewhat larger separ 
than the 8 pair. (1.8 A compared to 1.3 A 


2. Equipment 


Glass Flow System 


A glass flow system was employed to proy 
steady flow of sample vapor through the discl 
tube. A flow system decreases problems of cont 
nation by previous samples and simplifie s the pi 
ration of the sample for analysis. The 
shown photographically in figure 1 and schematix 
in figure 2, consists of the following components 

l. San ple bottles The sample bottles oe equip] 
with an inverted standard taper ground-glass jo 
so that rapid change of sample s can be made 

2 Sam ple connectors The 


syst 


connectors are stand- 


ard a ground-glass joints corresponding to those 


on the bottles. A special vacuum line (fig. 2 
vides a means of evacuating the air from the anal 
chamber before admitting sample into the main 
system 

3. Pressure regulator. In order to vary and cor 
trol pressure in the discharge tube a parallel networ! 


flow 


is provided, consisting of a capillary and a continu- 


ously variable metal leak valve. In most of th 
studies described in this paper, the variable 
leak valve was used because fine pressure 

and frequent adjustments were necessary 

ever, for routine measurements a fixed nonmetalli 
orifice is preferable. 

+. Pressure gage. ‘To measure pressure in the 
system, pressure line connects through stopcor Ks 
to various critical locations. Pressure measurement 
is made by means of a commercial ionization gag 
which utilizes alpha particles from a radium source: 
to provide ionization. The gage can be used from 
pressures of lu to lem. Pressures used in this report 
are gage readings. Absolute pressure of water vapor 
is found by multiplying by 1.163, the factor supplied 
by the gage manufacturer. 

Discharge tube. The discharge tube consists of 
glass inner cylinder surrounded by a glass wate! 
jacket. Water flows through the jacket at a constant 
rate controlled by a regulator at the tap. Air 
bubbles in the water, which might affect intensity 
measurements, are eliminated by an extractor shown 
at point G in figure 1. 

6. Pumps. The two main vacuum lines shown 
in figure 2 connect with independent pumping 
systems. ‘The left-hand line goes to a trap, single 
stage oil diffusion pump and fore ~pump while the 
right-hand line goes only to a trap and mechanical 
pump. The diffusion pump was used in 
investigations because a wide variation in conditions 
was to be studied. 


these 


metal 


control 


How- 








2.2. High-Frequency Source 


The energy necessary to excite the spectrum of 
the sample was provided by a commercial 150 Me 
150-w oscillator The energv is carried through a 
radio-frequency cable, which terminates in a tunable 
quarter-wavelength tank circuit with antennas held 
near the discharge tube. This terminal circuit is 
shown diagramatically in figure 2. A radio-frequency 
shield (point J in fig. 1) cuts down stray radiation 
which might affect other electronic equipment 


2.3. Monochromator 


The apparatus used to analyze the spectrum con 
sists of a high-resolution grating monochromator! 
with a photomultiplier detector and associated 
amplifier and recorde! \ detailed discussion of the 
original monochromator is given by Fastie [6] 

The plane diffraction grating used in the mono 
chromator is 3 by 3 inches and is ruled with 30,000 
lines per inch; the instrument has a prac tical resolv- 
ing power (A/AA) of 55,000 This resolution is 
sufficient to observe some of the fine structure split- 
ting of H, and D, lines. The choice of optimum 
resolving power is governed by two considerations 
it must be sufficient to resolve the lines to be analyzed 


This equipt t wa field t to the National Bureau 
ID 


Schematic diagram of flow system. 
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Height of H line (H,) is approximately 500 times that of H ghost line (G The 
zero line is arbitrarily taken 6.6 A from the H line, and intensities are measured as 
distances above the zero line. From left to right, traces are of 0.4-percent and 10- 
percent D)»O samples, respectively 


Frat Rb 3 Vethod of measuring intensities H " 


, 


and at the same time be low enough to allow for the 
integrated intensity of a line to be approximated by 
the peak height of that line. The resolving power 
was varied by varying the width of the entrance and 
exit slits until the 8 line of deuterium could be 
resolved in a normal water discharge (D/H=0.00015) 
(7). Figure 3 shows the H, D, and H ghost lines 
resolved, and indicates, on the of the half 
width of the line, a resolving power of 20,000. This 
resolving power proved suitable for subsequent 
studies. 

The first Rowland ghost of the 8 line of hydrogen 
was spaced about 2.0 A from the line and was ap- 
proximately one five-hundredths as intense. The 
8 line of deuterium appears between the H line and 
its first ghost. The ratio of ghost intensity to line 
intensity is a function of the optics of the system 
only, and under conditions of fixed optics may be 
taken as a constant. For low deuterium concentra- 
tions the ratio between the deuterium line and the H 
ghost is a fixed multiple of the ratio between the D 
and H lines. This fact proves useful in allowing 
comparison to be made between signals taken at 
the same scale of amplification, the calibration in 
this range being independent of the ratio of ampli- 
fication for different scales. 


basis 


2.4. Power Supply, Amplifier, and Recorder 


A chopper before the entrance slit of the mono- 
chromator converts the input signal to 60 eps. 
This signal is fed into a narrow band-pass syn- 
chronous amplifier, the output of which is recorded 
by a pen on 10-inch chart paper. The recorder 
has a response time of 1 second. The power supply 


for the photomultiplier is continuously variable f; 
500 to 1,300 v. By using the photomultiplier wit 
sensitive linear a-c amplifier, which has an accu 
amplification factor ranging from 2 to 800, \ 
small signals can be observed. 


3. Experimental Results 


The primary purpose of this study was to deter- 
mine the effects of numerous interrelated variables 
on the speed, accuracy, and sensitivity of the method 
The variables include pressure in the flow system 
flow rate, sample size, solute composition of the sam- 
ple, temperature, and oscillator characteristics 
These variables will be considered in detail. The 
percentage change in ratio per unit change in variable 
was computed for each operating curve and can be 
used to determine how closely each variable must be 
controlled to achieve a desired precision (see table 7 

The methods of measuring intensity ratios for the 
0- to 1-percent D,O range and the 1- to 99-percent 
D,O range are shown in figure 3. As previously dis- 
cussed, the deuterium to ghost ratio is used instead 
of the deuterium to hydrogen ratio in the 0- to 1- 
percent range. In the 1- to 99-percent range the D 
and H lines were measured directly at a lower amplifi- 
cation. In both concentration ranges the baseline, or 
background, is somewhat arbitrary and is taken at a 
fixed distance from the deuterium peak. This causes 
no error in the final measurement because the calibra- 
tion is made in the same manner as the measurement 


3.1. Contamination by a Previous Sample 


The speed with which a measurement can be taken 
is largely dependent on the time required to clean 
the system of one sample and saturate it with another 
so that a constant ratio is obtained. Metal parts 
are troublesome, but the metal leak valve usually 
was employed to gain fine pressure control in spite 
of its tendency to retain water with which it has 
previously been in contact. 

Figure 4 shows a series of experiments under ex- 
treme conditions in which a pure water sample was 
run and then switched to commercial heavy water 
Curve a shows the effect, under the slowest condi- 
tions, of using the metal leak valve with no interme- 
diate pumping or flushing. Curve 6 shows the effect 
of using a glass capillary tube and a glass stopcock 
for pressure control with no intermediate pumping 
or flushing. For routine measurement the capillary 
tube and stopcock were used under the following con- 
ditions. The system was pumped with and without 
sample between extreme pressures in the ratio of at 
least 15 to]. These changes were made rapidly for 
approximately 2 min, and then the pressure was 
stabilized to 400u. The effect of contamination by 
the previous sample was thereby reduced to that 
shown in curve ¢. 

These curves show that an all-glass system, with 
intermediate flushing, most quickly reduces contami- 
nation. Even with an extreme change in concentra- 
tion, equilibrium is reached to within 1 percent in 





30 100 300 


TIME , MINUTES 
Residual 
leak valve, with no intermediate pump 


glass stopcock, with no intermediate 
with intermediate pumping and 


Figure 4 contamination 


Pressure equals 400 uw. a, Meta 

g or flushing between samples 

imping or flushing glass stopeock 
Anab 


shing 
ss than 3 min and complete equilibrium in approxi- 
nately 1 hr. The flow system shown in figure 2 was 
designed for maximum variability, but for routine 
se, a much simpler system is advantageous. Also 
the elimination of side projections should reduce the 
time to remove the previous sample. 


3.2. Pressure 


The operating pressure in the discharge tube has 
an effect on the intensity ratio, especially in the 0- to 


i-percent D,O range (see fig. 5). Flow rate was not 
held constant, so actually the traces show ratio versus 
a combination of pressure and flow rate at constant 
pumping speed as determined by settings of the 
metal leak valve and the speed of the pump. The 
flow rate in milligrams per minute decreased from 
1.011 at a pressure of 900u to 0.0023 at 400, and 
0.001 at 200u. The pressure gage used had the dis- 
idvantage of large metal surfaces, resulting in appre- 
iable delay in removing traces of previous samples, 
in effect that is particularly undesirable when the 
oncentration is being changed frequently. There- 
fore, the gage was isolated from the system after 
nitial pressure setting and during the course of most 
ntensity-ratio measurements. The disadvantage of 
not taking continuous pressure readings was over- 
come by using the discharge length as a pressure 
ndicator because it is a sensitive function thereof 
\lthough this variation could not be used for absolute 
measurement of pressure, it provided a good check 
mn the constancy and reproducibility of pressure and 
flow conditions. 

Intensity as a function of pressure was investigated 
t various concentrations, and the results are shown 
n figure 6. The curves are quite similar at all con- 
entrations. An operating pressure of approximately 
{00u was chosen for routine use because it makes use 
f maximum intensity and the flattest portion of the 
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Figure 5 Intensity ratio versus pressure 
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intensity versus pressure curve. It also occurs at 
regions of the ratio versus pressure curves that are 
reasonably flat. The percentage change of ratio per 
unit change in pressure at 400u as computed from 
figure 5 is 0.020 percent per micron. This means 
that the pressure must be held to +5,y to achieve a 
precision of +0.1 percent. 
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3.3. Fractionation 


At 20° C there isa ams of 
between the vapor pressure of H,O and D.O, with 
HDO having an intermediate pressure [7]. Thus a 
sample that is being evaporated will begin to concen- 
trate deuterium To cher k on the error to be eX- 
pected from this effect, a 0.5-ml sample run 
continuously with frequent checks made 
on the D/H The results at an operating 
pressure of 200 uw are shown in figure 7. From this 
curve the minimum sample size can be computed for 
a given flow rate and allowabl This is shown 
in detail in section 5, 


2.3 mm of mercury 


was 
to dryness 
ratio 


error, 


3.4. Temperature 


and 
error 


The ratio of the vapor pressure of H,O, D,O, 
HDO changes with temperature, so that some 
might be expected from variation of sample tempera- 
ture. To study this effect, sample bottles were 
immersed in a constant-temperature bath during 
these determinations. Table | the results 
obtained at temperatures ranging from 0.4° to 25.8 
C. From this table the percentage change in ratio 
per unit change in sample temperature was com- 
puted to be 0.15 percent per deg C 


shows 


Tarsie | Intensity 1 { versus sample lempe 


3.5. Flow-Rate 


Flow rate affects the D/H ratio, the rate of frac- 
tionation, and the size of the smallest sample 
measurable. It is a function of pressure, the 
geometry of the flow system, and the speed of the 
pump. Table 2 gives data obtained for D,/H, ratio 


as a function of the flow rate at press 
This table shows that D,/H, « with 
creased flow rate and indicates that at least part 
the effect of change of pressure (fig. 5 
the reduced flow rate at lower 
percentage change of ratio per unit 
rate computed from table 2 is 0.5 
a minute. 


constant 


de creases 


IS Caused 
pressures 
change of 
perce nt per 


3.6. Oscillator and Discharge Characteristics 


The discharge characteristics can be varied i) 
changes in oscillator plate current, antenna tuning 
position of electrodes with respect to the tube. and 
the temperature of the cooling stream in the jacket 
of the discharge tube. Plate current is variablk 


between 100 and 200 ma, and the tuning is adjusted 
by the tuning capacitor in figure 2 as well as by th« 
position of the radio-frequency shield and electrodes 
The radio-fre quency shield and electrodes were plac ed 


to obtain maximum intensity, and the tuning capa 
tor was varied to give a range of discharge lengths 
Ratio of D,/H, as a function of discharge length is 
shown in table 3. The effect on D, H, of changing 
the oscillator plate current is shown in table 4. It 
that there is no appreciable effect of th: 
discharge length or plate current on the measured 
ratio. The temperature of the cooling stream was 
varied by changing the water flow rate through th. 
discharge-tube jacket. D,/H; varied than 
percent over the extreme range of temperature fro1 
a stream of the coldest water available from th 
tap (approximately 18° C) to no cooling (great 
than 200° C). 


is seen 


less 


TABLE 3 Intensity ratio versus percentage mazrimum 


le nath 
Pressure equals 500 » 
Percentage of 


maximum dis- 
charge length 


Sample 


100 
ss 
73. 
49 

100 
us 
91 
82.7 
75. 
4.7 


49 





3.7. Solute Composition 


\s the optical spectroscopic method of analysis 
extended for use with many aqueous solu- 
ys, preliminary investigation made of the 
riation of D, H, ratio with the amount of solids 
ssolved in the sample. This test was made with 
NaCl With concentrations of 0 percent, 1 percent, 
percent, and a saturated solution, no variation of 
D),/H,; was found within the experimental error. 


i\ be 
was 


4. Calibration Curves, Precision, and Repro- 
ducibility 


Calibration curves were necessary to achieve the 
chest accuracy Figure 8 shows the 0- to l-percent 
D.O range and figure 9 the 1- to 80-percent D,O 
* Samples were made by mixing portions of 
stilled water and commercial D,O (99.8%) that had 
been weighed to an accuracy of 0.1 mg. 

Table 5 gives examples of the isotope content ex- 
comparison, in ratio and in percent, 
ind shows the precision obtained under optimum 
conditions, 


inge 


essed, for 


Comparison of ratio method and percentage concen- 
i 


tration method of expressing isotope content 


Standard 
t euteri Coeftici 
. l 7 um deviation of . 
concentra ol 
average of 
ion Variation of 


t 
concentra 
100D/D+H tion © 


average * 


Ol 

ool 
ool 
Oooo! 
000007 
Qn0002 


jased on a set of 10 readis 


tandard deviation of average, expressed in percent 


bration curves were taken on samples up to 99.8 percent of DO (commer- 
y available), but no extensive measurements were made in the high-concen- 
mn range. The analysis in this range differs from that in the low ranges 
n the method of preparing standard samples. 
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2.0 


H,O (WEIGHED 


for 1 


range 


fion curve fa su 


FIGURE 9 


concentration 


to 9%-percent portion of the curve omitted because 


s were made uy 


21-percent 


Meas- 


was measured on a 


Reproducibility 
D,O sample and a 0.2-percent D,O sample 
urements were taken at varying time intervals, and 
in some cases with major changes taking place 
during the interval, both in the apparatus and in 


the concentration ratios Results are shown in 
figure 10 for 4 portions of the 21-percent sample. 








1 and reprod 


1).0 sample 


The data were taken over a period of 3 months, and 
dates are given from left to right 
on figure 10 


below, reading 


Sample 


Each point shows the average of 10 readings, 
with the precision depicted as vertical spread. The 
grand average is a solid horizontal line, and the 
standard deviation computed directly from all 21 
averages 1s shown by dotted lines on either side of 
the grand average Thus the relation between 
precision and reproducibility at this stage of develop- 
ment can be seen at a glance. The results obtained 
for the 0.2-percent sample were similar. Data for 


it are listed in table 6 


; 


Reproducihilit ; 


y of niens 


DO sam ple ove t l-mont 


April 28 
May | 
May 2 
May 4 
May 7 
Do 
Do 
May 8 
Do 


* (); is the measured intensity of the ghost 


f the H line 


5. Discussion 


The great number of variables that influence 
measurement of D/H ratios makes it impossibk 
be certain that the final operating conditions ar 
fact the optimum ones. Indeed, it is not necess 
to assume that there is a unique set of optim 
conditions for this analysis. The operating cor 
tions chosen therefore are the “‘ best set’’, based ul 
data available and the nature of the measurem: 
to be performed. Some of the final data used 
evaluate operating conditions are shown in tabl 
This table illustrates that there is no extr 
problem in keeping the operating conditions const 
enough so that the measurement can be made wit! 
the desired accuracy. 


TABLE 7 


Accessory factors enter into the choice of best 
conditions. The fractionation curve (fig. 7 
represented as RS* K, where R Is the isotope ratio 
S is the sample size, and k and K are experimental! 
constants In differential form the fractionatior 
curve can be represented as dP?/P k/S)dS, whe 
the left-hand term represents the percentage chang 
of ratio due to fractionation. Setting dR/R equal 
to e, the maximum tolerable fractionation erro! 
and realizing that dS fAt, the flow-rate times 
the time increment, we arrive at the following 
pression for minimum sample size 


Cail 


0.068 f At 


€ 


The value of &, 0.068, was obtained from the dat: 
shown in figure 7. As has been mentioned pre 

ously, flow-rate depends on operating pressure, flow- 
geometry, and pumping speed, and the time of 
measurement depends on the time required to cleat 
the system of the previous sample and to equalize 
the pressure. Hence, sample size is some complex 
function of the major operating variables. Under 
normal conditions at a 400-u operating press 

the flow rate was 0.0023 ml/min, and the time per 
determination was 20 min. The minimum samp 
size is therefore 0.00313/e, or 0.3 ml, for a 1-percent 
error. The absolute minimum sample size unde! 
present conditions involves an operating pressur 
of 200 « and a maximum error of 3 percent. Uniet 
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tion 


conditions, a 7-min measurement is possible, 
the mimimum sample size is 0.02 ml 
he operating conditions that were found to be 
suitable with the equipment described are 
yn in table 7. 
his study is not exhaustive, inasmuch as there 
many combinations of variables that were not 
d. In addition, no attempt was made to vary 
ematically such things as the oscillator frequency 
discharge-tube geometry. A few preliminary 
riments with a 2,240-Me oscillator indicate that 
able set of operating conditions can be set up at 
t frequency. 
\s pointed out in the beginning of this paper, the 
f object of these studies was to determine how 
minimize error, sample size, and time of determi- 
on. Of course, it is impossible to minimize all 
factors at Highest precision requires 
ny determinations per sample, hence a long time 
| large sample size. The minimum values attain- 
at the present stage of development and the 
ximum values obtained while minimizing the 
ers are Given as follows 


once 


Factor Minimum Maximum 


0. lto2 
cent 


0. O02 mil l 


per 3 percent 
mil 
' 60 mir 


With respect to precision and time per determina- 


the optical spectroscopic case compares favor- 


with other methods. It finds particular appli- 
tion where the water samples contain impurities 
Whereas other methods require one two 
purification, the sample vapor can be analyzed 
rectly by this method. This method appears to be 
msiderably more rapid and, for similar accuracy, 
ess expensive than other methods 
While it is not the function of this paper to discuss, 
detail, applications of hydrogen-deuterium analy- 
we shall review briefly in the following para- 
graphs some of the most important uses to which the 
optical spectroscopic method has been placed. The 
lilution technique for determining the total wate 
ontent of materials involved in these is de- 
ribed in reference Essentially, it consists in 
ssolving or suspending a known amount of the 
aterial to be analyzed in a mixture of D.O and 
HO and then determining spectroscopically the 
sulting change in the ratio of D,O to H,O. 
Determination of DO content of blood In this 
plication the blood, serum, or plasma is placed in 
sample bottle and the vapor analyzed in the 
ial manner. No prior preparation of the blood is 
essary and high precision is obtainable. This 
hnique finds application in the determination of 
tal body water by D,O dilution and in the study of 
iter exchange in capillary beds. Over 150 blood 
mples have been run, and no difficulty has been 
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or stages 


uses 


encountered in adapting the techniques to routine 
analysis 

2. Determination of water of hydration or water of 
crystallization The crystal to be analyzed 
dissolved in D,O, and the final intensity ratio gives 
a measure of the H,O originally in the crystal 
Several hydrates have been tried, and the results 
agree with other methods. 

3. Determination of labile hydrogen A weighed 
quantity of dry material containing labile hydro- 
gen dissolved in D,O and the resultant ratio 
measured. From the number of exchangeable hydro- 
gen atoms per molecule, information about structure 
of the molecule is obtained 

The most extensive work has been done with blood, 
for which a brief description of the routine applica- 
tion follows 

Army Medical Service gave injections of 100 g of 
commercially pure D,O to several hemorrhagic feve1 
patients and sampled blood at various time intervals 
The serum samples were analyzed at the National 
Bureau of Standards to determine the total body 
water content of the patients. Seventy-three sam- 
ples in the concentration range 0.029 to 0.425 percent 
vyave precisions of the order of 0 7-percent coefficient 
of variation of the average. Ten readings were 
taken per sample, using the metal leak valve for 
pressure control. The time per sample ranged 
between 25 and 55 min, depending on stability of 
operating conditions 

Reference [8] gives details of a study of 
capillary exchange rates of D,O and thiocyanate in 
the forearm of man. A study has since been made of 
the exchange rates from the artery of one arm to the 
vein of the other. In this case only five readings 
were taken per sample and metal were 
eliminated, thus cutting the total time 
down to 15 min, plus data calculation time 
precision of measurement was of the ordet 5 
percent coofficient of variation of the for 
17 samples in the concentration range 0.015 to 2.00 
percent of D.O. 


Is 


Is 


trans- 


parts 
per sample 

The 
of | 


average 
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Refractive Indices of Five Selected Optical Glasses* 
Robert E. Stephens and William S. Rodney 


indices of vVpical samples of five ont 
is were measured over a range of wavelet 


a WAY iat the values obtair 


ecise data on optical glasses should be readily & Son, Ltd., London Measurements in the infra- 
lable in considerable detail in order to satisfy | red and in the near ultraviolet regions were made on 
complex requirements of optical engineers and | each of two spectrometers. The first was of the 
gners. Only for the visual range of wavelengths | Wadsworth mirror type, improvised from a Zeiss 
refractive indices satisfactorily available for monochromator, the same instrument that was used 
er domestic or imported glasses, and almost no | for index determinations on magnesium oxide2 The 
are available on spectral transmittances second was the Gaertner spectrometer equipped with 
Ihe five classes conside red here were made and mnuirrors and a lead sulphide detector, the same appa- 
ealed in the Glass Section at the National Bureau | ratus that was used for index determinations on 
Standards under the direction of Clarence H potassium bromide 
hner, who has furnished the weight percentages The independently observed indices in the infrared 
wn in table 1, descriptive of their chemical | and also in the near ultraviolet regions were in fait 
position as computed from the batch composi- | agreement, and no systematic differences were de- 
With the exception of the rare-earth glass, | tected. The values resulting from the second set 
se are typical optical glasses and are probably Gaertner spectrometer should be somewhat more 
vy similar to the corresponding Bausch and Lomb | accurate because of a potentially higher accuracy in 
sses 7 1.517, v—64.5: n=—1.573, v—57.4: n angle measurement, and accordingly much more 
621, v=36.2; and 1.649, v=33.8, as listed in | weight is given to the Gaertner data in computing 
r catalogs under these D-line indices and v values averages 
Ihe refractive indices of these glasses were ac- The indices of relraction of these classes were 
ately determined in the visible region by means of | represented by dispersion equations of the form 
high-precision spectrometer made by E. R. Watts 


where X is in microns, and the parameters have appro- 
priate values for the various glasses as given in table 
2. The parameters were adjusted by the NBS Com- 


mpanion papers in this issue present values of the spectral trans 
properties, densities, and efficients of linear thermal expansion for 
fferent types of optical glasses produced by the Optical Glass Section of 
mal Bureau of Standard These measurements have been made in 
erent organizational parts of the Bureau appropriate to each type of 
; 


7 
ement, and are, therefore, presented separately For each type « 
have been taken from t me melt, and, therefore, these papers 


ollectively, a consister et lata for the five glasses studied 


glass 








TABLE 


omp ited sea of retrac dices ot glasse 8 


»F Rare eartt 
1.649 1.673 


8 “ 


1. 70407 
1. 7020: 


1. TO2R2 


1. 70250 
1. 60366 


HUS0. 


ARASO 
1. AS82272 


HS380) 


AS23686 
OS167 
1. 67038 


67618 
67400 
67386 


67322 
67027 


howe 


HORIS 
wo7n2 
HOT" 


6542 
68001 
65812 


65483 
6468 
65430 


65245 
65009 
‘4992 


H4839 
4508 
4128 
HAUS 





putation Laboratory, using the method of | 

squares. The magnitude and distribution of 

plus and minus residuals are such that the form 
are unquestionably adequate for representing 
data. The computed values, which are adjudge: 

be better than the observed, values are presente: 

table 3 

The necessity ot using a dispersion equation v 
six constants was not anticipated. Trials with { 
carefully adjusted constants, however, proved 
adequate in that for the first four glasses there y 
evident systematic differences between the obser 
and the computed indices. Such systematic fail 
for four-constant formulas were as large as 
10 10~° outside the visible region, and the deviat 
were more convincing in the infrared because of mor 
observed points there than in the violet. Moreoy 
the character of the systematic failures was similar 
for the different glasses, and thus the equation was 
considered inadequate in form. 

The addition of a term in \* was obviously a suit- 
able expedient for the infrared and one in 1/\ 
improved the ultraviolet end, and both seemed 
necessary. The resulting residuals are all less than 

5X10~°, with the exception of 2 (of opposite sign 
in the ultraviolet and, similarly, 2 in the infrare: 
One-half of all residuals are smaller than +0.7 10 
The remaining systematic trends between observed 


and computed values are not over +1 or 210 
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Spectral Transmissive Properties of Five Selected Optical 
Glasses* 


Harry J. Keegan, Marion A. Belknap, and Dorothy J. Cordrey ' 


Measurements of the spectral transmittance and internal transmittance of five optical 


glasses for the wavelength regior 


1. Introduction 


of a plate is defined as the 
o of radiant flux (P?.) transmitted by the plate 
the flux (?,) incident upon it. Internal trans- 
tance (7) is defined as the ratio of the flux (P 
ving at the second surface of the plate to that 
< (Po) penetrating the first surface. All values 
and 7, are functions of wavelength 
The transmittances and internal transmittances 
worted were determined for glasses of the NBS 
lts for which the chemical compositions are given 
Research Paper 2504 
lhe samples were prepared in the NBS Optical 
Shop in the form of cylinders or plates, 30 mm in 
imeter, with surfaces plane, parallel, and polished. 
The rare earth glass, n=1.673,? v=56.2, was pre- 
pared in thicknesses of 5 and 105 mm. The other 
ir glasses were prepared in a thickness of 1 mm 
addition to greater thicknesses. The thin samples 
ere compared against a blank beam, and spectral 
transmittance was so obtained. The thick samples 
were compared to the thin sample of the same kind 
of glass, the resulting ratios being taken as internal 
transmittance for a sample of thickness equal to the 
lifference between thick and thin specimens 


(he transmittance, 7 


2. Spectrophotometric Measurements 


Most spectrophotometers are designed primarily 
lor samples between 0- and 10- to 50-mm thickness, 
but can be used with samples of somewhat greater 
t The greater the thickness, however, the 
more the instrument is subject to errors resulting 
from refraction of the various rays in the beam, for 
although these rays are approximately parallel and 
at right angles to the surfaces, they are in general 
not strictly so. The result is, in effect, a change in 
P, to P, of the above defining relations for transmit- 
tance, P./P,, or for internal transmittance, P/P». 
(his is particularly true when, as in the present 


rhc kness. 


companion papers in this issue present values of the index of refraction 

ties, and coefficients of linear thermal expansion for five different types 

cal glasses produced by the Optical Glass Section of the National Bureau 

ndards These measurements have been made in the different organi- 
nal parts of the Bureau appropriate to each type of measurement, and 
therefore, presented separately For each type of glass the samples have 
taken from the same melt, and, therefore, these paners present, collectively, 

stent set of data for the five glasses studied 

Now Mrs. Frederick J. Does, E. L. duPont de Nemours & Co., Inc., Parlin, 
ndex of refraction for the sodium D lines, and », the Abbe value n—! 

. Where ng and nc are refractive indices for the hydrogen F and C lines, 
S was prepared in greater thickness of 97 mm and 
ared to the 5-mm sample. Values of 7; for 100 mm (for this glass are computed 
the values of 7;, for 92 mm so determined 


ectively 
ilass DF, n=1.049, » 


200 to 1,000 millimicrons are 


reported 


measurements with the thick samples, the compensa- 
tion beam of the spectrophotometer does not contain 
a sample of the same length (as it usually does, for 
example, in measuring the transmittance of a solution 
relative to its solvent 

Four types of spectrophotometers were available 
for this work. The designs of two of these instru- 
ments (the General Electric recording spectropho- 
tometer and the Kénig-Martens visual spectropho- 
tometer) are such as to give confidence in the results 
for thick samples. The other two instruments (the 
Beckman quartz photoelectric spectrophotometer 
and the Gibson photoelectric spectrophotometer [9})* 
were sensitive to beam disturbance with thick sam- 
ples. In the work on the rare-earth glass 1.673, it 
was found that the latter instrument was too sensi- 
tive to beam disturbance on insertion of the thick 
sample to give valid results, and it therefore was not 
used in the later measurements on the other glasses 

In the General Electric spectrophotometer (GE) 
(1, 2, 3] the beam of radiant energy, after transmis- 
sion through the sample, is incident upon a diffusing 
white surface of uniform reflectance forming part of 
the interior wall of a white-lined integrating sphere 
Any deviation or dispersion of the beam should not, 
therefore, cause appreciable error, if not so great that 
part of the beam impinges on the entrance aperture 
of the sphere. 

On the Kénig-Martens spectrophotometer (K-M) 
[4] the reverse but equivalent condition exists. The 
radiant energy is reflected from a diffusing surface of 
MgO, which is part of an enclosure completely cov- 
ered with MgO, except for the part containing the 
bases of the lamps from which the radiant energy 
originates. Multiple reflections should give a wall 
surface of essentially uniform radiance, and any dis- 
placement of this source in the photometric field by 
the thick sample should not introduce important 
error in the results 

The principal value of the K-M measurements was 
to confirm the reliability of the GE data even for the 
thickest samples. However, the precision of the 
visual measurements is less than the photoelectric, 
and except with the rare-earth glass 1.673, only a 
few check measurements were made on the K-M 
spectrophotometer. 

The Beckman quartz photoelectric spectrophotom- 
eter (model DU) (BQ) [5, 6] is somewhat susceptible 
to errors resulting from beam disturbance on inser- 
tion of the sample, and the thicker the sample the 
greater may be the error. On this instrument, the 


I 


‘ Figures 


lex references at the end of this 


n brackets indicate the inde paper 
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beam of radiant energy, after passage through the 
sample, is incident directly on the photocell. Al- 
though the window of the photocell is made diffusing, 
this is undoubtedly far from completely so, and any 
deviation or dispersion of the beam may well affect 
the response because of varying sensitivity of the 
photosensitive surface. The principal value of the 
Beckman spectrophotometer in the present work was 
to obtain the values of transmittance or internal 
transmittance in the ultraviolet, and it was not used 
at wavelengths longer than 460 millimicrons (mu 

Consistent with the above, the GE recording 
spectrophotometer was used over the range 400 to 
1,000 my, and the Beckman spectrophotometer over 
the range from 460 my toward shorter wavelengths 
into the ultraviolet as far as there was measurable 
transmission. The transmittance values obtained on 
the GE spectrophotometer were taken as a reference 
and the Beckman values adjusted to them when 
necessary in accordance with the following pro- 
cedure 1) The GE transmittance values were used 
from 460 to 1,000 mau, subject to an occasional 
smoothing not exceeding 0.001; (2) from 410 to 460 
mu the GE data were used to compute an adjust- 
ment factor to apply to all of the Beckman data 
This adjustment factor was the average of the ratios 
(Tor/ Taq) at each of the six wavelengths from 410 to 
160°'my. In this range the adopted values were then en — 
based on both the GE data and the adjusted Beck- WAVELENGTH, mp 
man data; (3) from 400 my toward shorter wave- Figure | Spectral transmittance of 
lengths the adopted values were based on these glasses made at NBS for the ultra 
adjusted Beckman data Sree eae 

All the values by the Beckman spectrophotometer ° ’ i 
were obtained with the samples at 25°C. Values by 
the GE spectrophotometer (visible and infrared 
were obtained at a room temperature close to 25° C 


>= 


and should be valid also at 25° C 


TRANSMIT TANCE 
Imm 
n= 517 
L574 
L620 
L649 


n 
n 
n 


TRANSMIT TANCE , % 


3. Spectrophotometric Results 
INTERNAL TRANSMITTANCE 
Curves of the transmittances for 1-mm thickness a. 
for four glasses are shown in figure 1 (numbered 1 to 
4), together with the transmittance curve for the 
5-mm thickness of the rare-earth glass 1.673, num- 
bered 5 in the figure 

Curves for the internal transmittances for 4- and 
100-mm thicknesses of the same four glasses are 
shown in figure 2, together with the internal trans- 
mittance curve for the 100-mm thickness of the rare- 
earth glass 1.673. 

The values of transmittance and internal trans- 
mittance are given in table 1. All are subject to 
slight uncertainty in the last figure. Examination of 
the values of 7; for 100 mm for the rare-earth glass | 
1.673 shows traces of the absorption bands near 580 
and 740 mg commonly found in didymium glasses. 


TRANSMITTANCE , % 


INTERNAL 


600 . 
WAVELENGTH ,m 

* This indicated procedure was followed for 4 of the 5 glasses. For rare-earth 4 : 

Oye 1.8, a slightly cea pF nye nae was used Fraure 2. Spectral internal transmittance of five 9% 
t the steep portions of the curves (figs. 1 and 2), the transmittance and y . ¢ , a } 
S for sibl 

internal transmittance vary importantly with temperature, in these cases in the opti al glasse s made at NB Jor the ultra tole t, waite, 
ultraviolet region of the spectrum [7} near-t nfrared spectrum 
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Values of 7’, for 4 mm are not tabulated above 400 4. References 
or 450 mu where the values exceed 0.995 As there ; ' H ; — 
is uncertainty of 0.001 or 0.002 in such values, it | |!) 4 ‘ lardy istory of the design of the reg 
, . spectrophotometer J Opt. Soc Am. 28, 360 (193 
would be preferable, where the absorption Is sO 2) J. L. Michaelson, Construction of the General | 
slight, to compute the values of 7 for 4 mm (or for recording spectrophotometer J Ont Soc Ay 
| orl I I | 
any other small thickness) from the values of 7, for 365 (1938 
100 mm than to measure them directly 3] K. 8. Gibson and H. J. Keegan, Calibration and oper 
>. o . of the General Electric recording spectrophot« 
lhe slight error in 7; resulting from the multiple of the National Bureau of Standards, J. Opt. So 
reflections through different thicknesses of glass [8] 28, 372 (1938 
is ignored in these results. For either high or low | [4] H. J. MeNicholas, Equipment for routing 
transmittance the error is of the order of 0.001 transmission _ and reflection measurements, 
. Research 1, 793 (1928) RP30. 
In view of the equal or greatel uncertainties indicated 5] H. H. Cary and A. O. Beckman, A quartz photos 
above computation of such corrections was con- spectrophotometer, J. Opt. Soc. Am. 31, 682 (194 
S. Gibson and M. M. Balcom, Transmission mea 
ments on the Jeckman quarts spectrophotomet« 
Research NBS 38, 601 (1947 tPI798: also J 
Soc. Am. 37, 593 (1947 
K S. Gibson, Spectrophotometry NBS Circular 
ial . ; 1949 See figures 29 and 30 
rhe authors express their appreciation to K. S. | (9) 7. Smith. Note on measurements of glass absort 
Gibson, Chief, NBS Photometry and Colorimetry Proc. Phys. Soc. 58, 473 (1946 
Section, who coordinated the work of these spectro- K. S. Gibson, Direct-reading photoelectric measure 
, of spectrs ‘ans ssion J ) oc m. and | 
photometric measurements and results; to Kenneth : , ee pt. Soc. Am. and 
— > ee ns a, OFTO a0 
L. Kelly for his assistance in the computations; and 
to John C. Schleter for the preparation of the graphs WaAsHINGTON, November 19, 1953. 
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Densities of Five Selected Optical Glasses* 
Charles T. Collett 


five optical glasses were ds 


ted by type of glass, index of refraction, and 


measurements 
n the form of a cylindrical disk 30 mm in dia- 
r and 5 mm thick. The observed densities of 
samples were derived from their masses, deter- 
weighing in air, and their volumes, deter- 
hydrostatic we ighing, with corrections for 


glass used in these 


ich sample ol 


ed by 
d by 
vVanecy 
or the hydrostatic 
led from the left pan of the 
passed down through hole in the case into a 
net insulated against changes in temperature 
air currents. The immersing liquid was dis- 
d water in a Dewar flask closed by a stopper, 
igh which holes of small diameter admitted the 
ension wire, a glass stirring rod, and a platinum 
stance thermomete! \ calvanometer and Muel- 
type resistance bridge continuously indicated the 
stance of the thermometer, from which the water 
perature was calculated. 
Che effect of surface tension was reduced as much 
by the use of gold-plated platinum wire 
mm) at the region of inter- 
on with the water surface. The errors caused 
surface tension were further reduced by 
on between the two steps of the weighing process 


was SUS- 


weighing, a basket 


balance by a wire 


ossible 
small diameter (0.20 


cancel- 


described 
\ constant load remained on the right pan of the 
ince, and the immersed basket was at all times 
pended from the left pan. The first step was to 
ce each glass on the left pan with weights, W,, 
ficient to produce equilibrium. For the second 
p, the sample was transferred to the basket, and 
ilibrium was restored by increasing the weights 
the pan to the value 
Provided that W; and W, 
change of density occurred 
volume of the glass sample, V,, 


were determined while 
in the water or air, 
was given by the 


ation 


(W.—W (i—es ‘ 


Pr Pa 


V; 


npanion papers in this issue present values of the index o 
tran smissive properties, and coefficients of linear thermal expans 
different types of optical glasses produced by the Optical Glass Section 

National Bureau of Standards measurements have been made 
lifferent organizational! parts of the 
irement, and are, therefore, presented separately For each type 
samples have been taken from the same melt, and, therefore, 


collectively, a consistent set of data for the five glasses studied 


These 
Bureau appropriate 


termine 


{ refractic 


to each type of 
of glass 
these papers 


d by the hvdrostatic we 
Abbe value 


where p and Po represents the densities in grams per 
milliliter of air and water and 8.4 the density of = 
brass weights. The air density was de termined | 
means of a special buoyancy balance 

Weight@iwere observed to 0.1 mg, and temperatures 
to 0.001 deg ( The observed values of density are 
shown in table | ; 


Winkelmann and Schott, and Russ,’ among others, 
made density determinations on selected glasses and 
derived factors for various component oxides By 
use of these factors, the density of the glass could | 
roughly computed 

The last two columns in table 1 give density values 
thus computed The differences between these and 
the observed values vary from less than 0.5 percent 
than 7 percent. This is not surprising, 
in view of Hovestadt’s * comparisons of 29 glasses, 
where trial density computations based on 
factors revealed differences as large as 10 percent. 


to more 


these 


N. 8. Osborne ely nd earce, Bul. BS 9,32 


see Pp. 378 

2 Geo, W. Morey, The pr 
New York, N. Y., 1938 

H. Hovestadt, 
(Macmillan & Co., Ltd., 


nhold Publ 


Jena glass and its scic 
London, 1902 


WasnHinoton, August 14, 1953. 











Thermal Expansion of Five Selected Optical Glasses 


Peter Hidnert 


Coefficients of linear thermal expansion of five optical glasses were determined by a 


strall 


gare 


of expansion of two crown glasses are 7.6 10~* per deg ¢ 
but the coefficients of expansion of two flint 
larger than the other glasses 


have a slightly lower coefficient of expansion 


glasses are 10 percent 


his paper gives coefficients of linear thermal 
insion of five samples of optical glasses. The 
nical composition of these glasses are given in 
arch Paper 2504 
he diameter of each sample of glass was 30 mm 
length of each sample was 105 mm, except 
ple DF (table 1), which was 97 mm long. Deter- 
tions were made on the linear thermal expansion 
hese glasses on heating and cooling between 
and 60° C An electric-resistance strain-gage 
hod was used. When a strain gage, consisting 
fine wire, is attached to a sample of solid material, 
nges in length of the wire accompanying thermal 
ansion of the sample proportional to the 
nges in electrical resistance of the wire. Infor- 
on about commercially available types of wire 
n gages has been given by Ramberg.’ 
(wo SR4 (type A 11) strain gages connected in 
es, were attached with a cellulose nitrate cement 
opposite sides of each sample of glass and of 
h of three reference samples. These reference 
ples were duplicates of samples of invar, glass, 
steel on which coefficients of linear thermal 
insion had previously been determined by the 
ision micrometric method* of the Bureau. 
er the cement was sufficiently dried, the cemented 
ves were coated with wax. The glass and refer- 
samples with their attached gages were then 
bjected to two cycles of heating and cooling 
tween room temperature and 60° C, in order to 
low the gages to adjust themselves to the samples 
before observations were made. During these cycles 
f heating and cooling, all samples were held at the 
ximum temperature for 2 hours. 
The samples of glass and the reference samples 
vith attached strain gages were supported vertically 
an electrically controlled temperature cabinet. 
Two copper leads extended from each sample and 
from each reference sample to an SR4 strain indi- 


are 


this issue present values of the index of refraction, 
and densities for five different types of optical 
s produced by the Optical Glass Section of the National Bureau of Stand- 
These measurements have been made in the different organizational parts 
Bureau appropriate to each type of measurement, and are, therefore, 
ted separately. For each type of glass the samples have been taken from 
sme melt, and, therefore, papers present, collectively, a consistent 
{ data for the five glasses studied 
Ramberg, Contributions of electricity to mechanics, J 
Sept. 1949 
Hidnert and W 


¢ companion papers in 
i] transmissive properties 


these 
Wash. Acad. Sci 
38, 281 


Souder, Thermal expansion of solids, NBS Circular 484 


294330 —54 


method for the range between room temperature and 60° (¢ 


The coefficients 


A rare-earth glass was found to 


cator (type K Four thermo- 
couples placed symmetrically in the region of the 
glass samples and the reference samples were con- 
nected to a precision potentiometer 

Observations were made at room temperature and 
at 60° C, during three cycles of heating and cooling, 
on each sample of glass in conjunction with each 
reference sample of known expansion. These obser 
vations at constant temperatures, consisted of strain- 
indicator readings recorded to the nearest microinch 
per inch and of potentiometer readings recorded to 
the nearest microvolt. The latter readings were 
converted to temperatures to the nearest 0.1 deg C 
The coefficients of linear thermal expansion of the 
five samples of glass given in table 1, were derived 
from the data. These coefficients of expansion are 
believed to be accurate to +0.1* 107°. 


copper-constantan 


linear thermal expansion of 


( al glasse s 


TABLE 1 Coe fli of 


opt 


¢ No factor was give 

The last column in table 1 gives coefficients of 
expansion computed from English and Turner's * 
factors.’ One-third of the sum of the products of the 
factors and the precentages of the corresponding 
constituent oxides of a glass gives a computed coeffi- 
cient of linear thermal expansion. The agreement 
between the observed and computed coefficients of 
linear expansion of the glasses is good. 

‘S. English and W. E. 8. Turner, Relationship between chemical composi 
tion and the thermal expansion of glasses, J. Am. Ceram. Soc. 10, 551 (1927 
Corrections on “Relationship between chemical composition and the thermal 


expansion of glasses,"’ J. Am. Ceram. Soc, 12, 760 (1929 
' For computing coefficients of cubical expansion 
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The author acknowledges his indebtedness to 
Valore, Jr. of the NBS Structural Engineering S, 
for helpful suggestions in the application ; 
strain-gage method and to H.S. Krider of the Ly 


The coefficients of expansion of the two crown 
glasses (BSC 1.517 and LBC 1.574) and the rare 
earth glass 1.673 are practically the same and about 
10 percent smaller than the coefficients of expansion 
of the two flint glasses (F 


1.620 and DF 1.649). | Section for assistance in the experimental work 


WasHINGTON, July 14, 1953. 
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E valuation of the Exponential Integral for Large Complex 
Arguments 


John Todd 


Two methods of evaluating the exponential integral for (large) complex arguments are 
discussed. It is shown that the Laguerre quadrature method is more efficient than the 
asymptotic expansion. Examples are given to show the practicability of the Laguerre 
method, 


1. Introduction 


An extensive table of the exponential integral 


E\(z) | (e~“/u)du zs=a2+uy 


has been prepared by the National Bureau of Standards [1];' the introduction to the table 
gives a precise definition of this function. This table covers the range |z| <20, |y| <20, with 
arguments variously spaced. In order to compute /\(z) outside this range, (or within this 
range, at points where interpolation is awkward, which are those with comparatively large 
argument), several methods can be considered. We examine two in detail and recommend the 
Laguerre quadrature, which is certainly practicable in the case of isolated arguments 


2. Asymptotic Series 


The first of these is the use of the asymptotic series (see e. g., [1] 


E\(z)~« ‘(--5 


The remainder after n terms in the expansion for e’/(2) can be obtained by integration by parts 
and is 


If we put u 


which can be estimated as 


if <0, 
from which the asymptotic character of the series is evident. In fact, 
n+1 
, , 
(z*+ y*)? 


and so S, decreases as n increases until n exceeds (z*+y’)'. To obtain the least the least value 
of S, we have to take n==(z’+y’)'=|2|. For this value of n, we find, using Stirling’s formula, 


n 


In the second case, the least value of 7, is formally the same, but 


that S, is about (22/n)* 


| Figures in brackets indicate the literature references at the end of this paper. 
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now n—-ly|. The decrease in the effectiveness of this representation as z approaches the 


negative real axis is apparent from the estimates 
We note that the evaluation of the sums of the first n terms of the series for e*F, (2). as a 
polynomial in '. by successive multiplication by 2~' for instance, involves about 4n real 


multiplications, z being complex 


3. Laguerre Quadrature 


The second method is the use of a Laguerre quadrature. We shall show that this is mor 
efficient than the asymptotic series method just described 


We have 


where 


It is well known [3] that, for any n, 


. 


/ | e~*. f(t)-dt=-Sr\ f(z™)=Q, 
0 


where the z{"’ are the zeros of the Laguerre polynomial L,(t), and the \{*’ are the correspond- 
ing Christoffel numbers. The x’, have been tabulated for n=1(1)15, i=1(1)n by 
Salzer and Zucker [2] 

Moreover, it is known [3] that for some £, 0<&< @, 


t=|J—Q] 


We shall now estimate this error. It is convenient to handle the two integrals separately. 


We calculate the 2nth derivatives of 


> ¥ ol ° V] 
f(t) " 7) f(t) ~- 
(7 +f)*-+ y° z (7 +-t) 


and find 
(n!)*r7** -' cos (2n+1)6,, 


(n!)*r>7* -* sin (2n+ 1)62, 


where the £,—£,(z, y) satisfy 0<£ *-+-y*|', and 6, is defined by cos @, 


sin 0,=y/r 


The difficulties that were noticed in section 2, as 2 approaches the negative real axis, are stil] 
present. As before, there is an optimum value of n, which is about || in the first case, |y 
in the second. With this choice of n, the value of our bounds for the R, are about 2xe-2" in 
each case. 

The important fact to notice is that our estimate for the error committed by using an 
n-point Laguerre quadrature is approximately the square of that for the estimate of the error 
committed by neglecting all terms in the asymptotic series after the nth. This, together with 
the fact that about the same number (4n) of real multiplications are required in both methods, 
indicates that the Laguerre method is more efficient. 
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4. Examples 


In order to demonstrate the practicability of the Laguerre method, and to show that our 
conclusions are not dependent on the crudity of our appraisals of the remainders, we give three 
examples in detail. We choose n=5, which is optimum in one but not in the others. For 
completeness, we repeat (from [2]) the values of z,, \,“ 


= Ya 
. 26356 03197 18 0.52175 56105 83 
. 41340 30591 07 39866 68110 83 
59642 57710 41 07594 24496 817 
7.08581 00058 59 . 00361 17586 7992 
12. 64080 08442 76 . 00002 33699 723858 


The superscripts (5) will be dropped from now on. 
We shall work to about 10 decimals throughout our examples. We begin with a case of 
real: z=10 and follow with z=1+ 10i, and then z 10+-57; it is only in the last case that 
we are, with n=5, near the optimum number of terms. 


10. 
In this case f,—0, and we have only to deal with 


I, : ee ce A TAG: Q, 


“(10+0) —_ 10+2, 


A complete worksheet for this calculation follows. The sum Q; is cumulated on the machine, 
there being no reason to record individual products 


a,=— 10+ 27; az" 
26356 03197 0.09743 20771 
. 41340 30591 . 08761 62872 
3. 59642 57710 . 07354 87412 
. 08581 00059 05852 81002 
2.64080 08443 . 04416 80490 


0, > Aa: '—(), 09156 33319 


The correct value of e'®£,(10) is 0.09156 3334. The actual error is about 210 


whereas our estimate is 144X10-*. Using five terms of the asymptotic series, we obtain the 


value 0.09164. There is thus an actual error of about 810~°, compared with an estimated 
error of 12X10 
b 2 l t 10% 
In this case we have 
et F(1+-100)=—h-th, 


where 


-dt=S*) 


Aq’ 
— (1 


Lz,?+1008° 








A complete work sheet for this calculation follows. It is convenient to compute the common 
factors y,, in braces above, and then obtain the sums Savy, and 5)10y, by accumulation on 
the machine 

a l+2, 8,=100+- aj Vi=Ai/ Bi 

26356 0320 101. 59658 47 0.00513 55625 

41340 3059 105. 82451 43 .00376 72444 

59642 5771 121. 12712 99 . 00062 69648 

O8581 0006 165. 38032 35 . 00002 18391 

64080 0844 286. 07144 76 00000 00817 

Re FL + 101) SS avy; = 0.01864 0471 


E\(1+-10i)= — 35107, = — 0.09551 6925 


The correct value of e'*'F,(14+-10i) is 0.01864 049—0.09551 6887. 
The actual error is of the order of 5 107°, whereas our estimate is about 14410-*. Using 


five terms of the asymptotic series gives 
0.018595 — 0.0956002. 


There is thus an actual error of about 1 10-*, compared with an estimated error of about 
1.2 107° 


a lO+z B,= 25+ aj ¥i=Ai/B; 

9. 73643 9680 119. 79825 76 0.00435 52855 
8. 58659 6941 98. 72964 70 00403 79645 
6.40357 4229 00576 29 .00115 05428 
2.91418 9994 33. 49250 33 . 00010 78378 
2.64080 0844 31. 97382 91 00000 07309 


-10+ 5i)= S35 avy: = — 0.08475 7264 


‘Ey (—10+-51)=— 335 0.04826 1807 


> 5 
The correct value of « os 10+-52) 1s 


08475 749—0.04826 039%. 
The actual error is of the order of 2 10~-*, whereas our estimate is about 310-*. Using 
five terms of the asymptotic series, we obtain the value 


0.084855 —0.0482771. 


There is thus an actual error of about 110~*, compared with an estimated error of about 
8x 10-° 


5. Remarks 


1. A third method of attack is conceivable: The integration in the complex plane of the 
differential equation for £,(z), taking 2~' as the independent variable, extending the ideas of 
Fox and Miller [4]. 

2. Itis conceivable that the Laguerre method would be effective in cases when the integrand 
is not precisely of the form e~‘f(t). A case of importance might be that of the generalized 
exponential integral ‘. (e~”/w)du, where w=(a?+u*)'. In such a case it would be conven- 


ient to follow the device of A. Reiz [5] and use the quadrature in the form 


f, Poa Dur Fe), 
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where »{ =!" exp z{; the »{” have been tabulated by Salzer and Zucker [2] 


3. The table of Mashiko [6] gives /,(z), or rather auxiliary functions from which this can 
» obtained, for re’ when 7 0.01(0.01)0.2 and @=0(2°)60°(1°)90 
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On Nearly Triangular Matrices’ 


A. M. Ostrowski * 


is presented of the change in the inverse of a triangular matrix if on one 
the other side the non 


replaced by sufficiently small numbers and o 
sufficiently small amounts 


s Ss1on 


zeros are 
are varied by 


gy elements & 
1. Introduction 


Consider a system of linear equations 


- n) are not equal to 0 and 


Me 
\/ the inequalities 


A, in which all diagonal elements a,, 


with the matrix 
the elements off the diagonal satisfy for two positive numbers m 
ay»| S mia,, 


a,,| = Mia,, 
triangular’ 


If m is very small, the system does not essentially differ from the corresponding 
uw are replaced by zeros and the matrix of which will be denoted 


system in which all @,, with »- 
It then appears plausible that the solution of this triangular system does not differ 


by A 
very much from that of the system (1 


However, the value of the determinant of the order 7 


10 


shows that if ./ is, for instance, greater than or equal to 10, the determinant of our system 
\ detailed study of the problems 
appears, therefore, to be of importance and 


will not be necessarily different from zero unless m 


connected with the matrices characterized by 


interest 
As the first problem in this connection, we give a necessary and sufficient condition that any 
If m< M, this condition is given by 


matrix A satisfying conditions (2) be nonsincular 


M 
M)"’ 


m 


m 


This paper was prepared under a National Bureau of Standards contract with American Ur 
American Ur rsity and University of Basle, Switzerland 
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and, if m V/, by 


In order to obtain a precise measure of the influence of the change from A to A we have to 
discuss the estimates (for convenient norms) of the norm of the matrix A A 

We consider, in particular, two such norms defined in section 5 and denoted by 
A A » & which are particularly suitable for the problems of numerical analysis 
Assuming, without loss of generality, that @,,—1 («—1, .n), we show that for given 
values of m, Mf such that 72 1.5/n, n2=4, we have 


V/ 


where 1—é6 is the smallest modulus of the determinant attainable for the matrices A and is 


connected with m by the relation 


Vm 


mm 
4 


V/ 


V/ , nM 


M, Vl 


If MV i/n, the formula (4) need not be valid any longer, but we can prove in this case 
the relation 
onm 


nm 


is valid as long as m remains less than 1/2n. 


The estimate (7) is not a “‘best”’ estimate for all values of W<1.5/n, but still it is not far 
from the best, since for m=.\/<1/(n—1) we have 


which cannot be improved for any value of m< 1/(n—1) 


The condition (3) is derived in section 4, theorem I However, we derive it as a special 


case of a more general theorem, where in the inequalities (2) the expressions m, Mf depend on u, 


that is to say, change from one row to another. The necessary and sufficient condition for all 
matrices A to be regular (theorem A, section 4) is in this case rather unwieldy, but still may be 
very useful in some cases because it contains 2n—2 instead of two essential parameters. The 
direct derivation of theorem BP is, of course, much simpler, since, as the reader will immediately 
see, the computations of the determinant @, in section 2 can be considerably shortened in this 
case. The connection between the formal algebra of sections 2 and 3 and theorems A and B 
is provided by a result concerning the so-called /7/-determinants and .M-determinants published 
16 years ago [1] The results about the norms |A~'—A~'|, are obtained by using the 
explicit representation of the inverse matrix of a certain matrix A,, which provides a majorant 
for all matrices A~'. The formulas giving 4>' are derived in the second part of section 3, and 
in section 6 the norms A;'—A°~'|, are derived and discussed. The corresponding inequalities 
for |A~'—A® >", are then obtained in section 8, by using a new theorem (lemma IIT) concern- 
ing the connection between the //-determinants and the M-determinants 


> This condition (3¢) is alread: yntained in some results in a paper [2]‘, and also can be deduced from a well-known theorem of the theory « 
leterminants discussed in [4 


* Figures in brackets indicate the literature references at the end of this paper 
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In this section is given another application of this theorem in estimating the variation in 
the inverse matrix of a triangular matrix satisfying the conditions (2) with m=0. We obtain 
an unexpectedly simple and elegant formula (120 


In section 8, the results are applied explicitly to the problems concerning the linear system 


It may be finally remarked that these results, with obvious changes, remain valid if in 


the matrix of (1) the rows and the columns are interchanged. although no mention is made of 
it explicitly at every step 


2. Value of the Determinant °, 


Let K bye ce fined by 


where in the u-th row all elements to the left of the main diagonal and 
equal to x,_, the next element to the right is 6 
are here independent variables 


on this diagonal are 


and all other elements are 0. «,, 4, 
In subtracting the second column from the first, we 
kK K,»—6,)K and therefore the following formula, valid also for n= 1, 2 


obtain 


K, 


Consider now for n >3 the determinant 





We have in particular 


Developing 7, in the elements of the first line and using the value (10) of K,,, we obtain 


> 4 


Bw Y. 


nu (1 4 


r B, Y2 8 
Is )-3 


6.6 0 
we obtain 


~ Bn" Y: 
> . un ( 2 ), 


where II is identically 1, and therefore finally 


If we now put 


this becomes 7), if the indices of 8,, y,, 6, are replaced by their complements with respect to 


n-+-1. We obtain then from (14) 








T* 
and finally, if in both products v is replaced by n+ 1 
T?=—-> 8,16 0 


Consider now for n=>3 the determinant 


m m m a 
If we subtract here the last column from each of the preceding ones, we obtain 
V/ 0) 
Vi, iy \f,4 
Vv M, 


V, —_ * ' M, 


aA» MN, a, h,!) (\Q, Tt 


Here the subdeterminant corresponding to ’ element of the last row is obviously 


Il (.\4,+.a,), so that Q, is the sum of a, iT (\f,+a@,), and 


M+ 
M, 
M 











This last determinant becomes 7,,* if we put 


é, V/, a, 7, M, m,, 


and has therefore by (16) the value 
SM,’ 
>> M, 


We obtain therefore for 2, the expression 


2,=a, ft (M.+a)—3°M, I (M4 
- l AN I 


« 


The determinant @, can be also written in the form 


Vv, 


Vi, 


and we obtain therefore from (18 


n 
Q, =a, I (m,+ a,) 


rar 


I | T Oly nt (.M, + at,) 


3. The Matrix A,, and Its Inverse 


We consider now the matrix 


M M 


M 





x m m 





m 


Its determinant is obtained from Q, in putting in (20) 


a) =a l; Mo=M™s, m=: 
We obtain 


n 
A,| = (m+1)*"'—m>*\(m+1)"-*(M+1)* 





and this becomes, if VM #m, 


Li M(m-4 


°s|“ Fi — am 


\/ we obtain from (22 


while for m in letting 


A, |} (n ] mi (1 


In particular, if 0<m- 


m 


We assume now in particular 


If we introduce the abbreviations 


M,= (ite 


we can write (22) in the form 


M—wm)i\A, M(i+nm+mm,)—m(1 


and therefore, if we put 


6 M, 
M 


m M 


It follows from (25) for n 


h M: == 
( Vv ) M—m 


we have 


U. Mm» a 
M 


ae 
m y=? 


so that from (28 


M, 


m 


M 


It follows in particular that if (24) holds, then 6 


In solving (29) with respect to m and 4, we 
6 


= 
Mit yy 


0 
M, 


M » ti 
m 


M 


M,m < 6 
| - 


also obtained from the formulas given in the proof 


t nM 


V/ 


VW). 


+m)" (7 


V/, a necessary and sufficient condition for |A, 


M 
i n 


MM, Vf mM(M, 


m 


A, |, 


nm, 


m 


y(t) a 


obtain 


VM 


of theorem III in [3a, p. 113) 





Weal. 











In order to discuss the meaning of the condition (24), consider the curve 


For « >0 the curve (31) has one maximum at z= 1/(n—1) and an inflection point at «=2/(n—1) 


In figure 1 the curve (31) is drawn (computed for n=5 The portion of this curve from 


the point 0 to the highest poml T will be denoted as the asce nding branch and the portion 
between 7 and 4 2 as the descending branch 


f(x) 








| I 
1 ! 
| | 
M2 
n-! 
x 
f(x) = (1+x)> 


Fiaure | 


If we assume that in (24) m is less than MM, we see that either 


In order to find for a given \f>1/(n—1) the range of values of m satisfying (24) we find 
on the curve the second point Py with the same ordinate as the point s=.//; then if the abcissa 
m(M) corresponds to Py, we have 0<m< m(M) 


Finally we prove from (24) and m< M that it always follows that 


yn 
mi" 


Indeed, in proving (33) we can obviously assume that .V is greater than | and therefore 
the descending branch in the diagram. But then we have 


Aq*-"' M 
) 1+ M*-! + Af’ 
} ” aati 


1+ 


and therefore in raising it to the nth power 


M 
Y “(+M)" 


1 


We see that 1/\/"~' must lie, in figure 1, between Mf and m(.M), and it follows that 1//""'>m, 
that is (33) 





We are now going to determine the inverse matrix 
A ( Qu») (uv 


of A,, if its determinant (22) resp. (23) is not equal to 0 


the corresponding linear system 
MW » > zr 
—_ a 


for indeterminate 





Then the system (35) is equivalent with 


l Ma, 


or in putting 


From (37) for uw=n we have bv (36) 


‘ ’ 


ay | 
y= 


where for m V/, q 1 the coefficient of Xs is to be replaced bv n Since for m#M 


l m M q’ Mq a m 
l m M M mn ; 


we have from (39) and (22) 


Q >> aq" ’ mM), 


(mn = M), 


291330 54 











while for m= .M we obtain again (40) in defining Q by 


m-+-1 (m-—+-1)’ 


7 
l \vi l)m 4. 


If we now replace in y, by ¥s +a, and use for s the expression from (40). 
| 


l S a,q” 


But here the expression within the brackets is q*~', and we obtain therefore 


‘ 


a ag” +a, AV aq" , 


am, Y 


and finally in introducing a, from (36 


l — AV aq" m 
\ ey s 


1’ M4+1° Me+1 


Introduce the two following triangular matrices 


. 0 














0 

(? 

Loe : (u<_pv) 
We see then that in (42) the matrix corresponding to the first right-hand sum is 


q l I m M 
M+1~" (M+1) 
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whereas to the second right-hand sum in (42) corresponds the matrix 


Qe. 
M+it ( 


where -, is the unit matrix We obtain 


m M MQq" \u, 


rl l MQq ; MQq" 
(M+1)?* (M+1» ' )E.+ 


A 4 
Cy | M-+1 M +1) 


The coefficient of U, is here equal to (m@)/(\f+-1)*, both in the case (41) and (41°), and we 


have finally 
mq l WUQ4¢q VE VWQq' 


wie Haat (iy)! tari a?) 


We now denote by s, the sum of the elements in the uth row of the matrix (47) and by ¢, 
the corresponding sum in the uth column. If we first assume m#.M/, for the matrices U 


and V’,, the wth row sums are, respectively, equal to 


and the vth column sums are 


We obtain now from (47) 


VWQq’ WqQ ! 
4 


( 
€ Imag” '—_m-+ Mq’ Mq*~' Vq VUq*'|\+ M+1 


M +1 


m= M @ |. Mq’ m V/ q” 


and this is by (41) equal to (M/+-1)(Qq* We have finally 


and this remains true also for m=, as is immediately seen 


Similarly, we have 


(VW t ] *t, m@ T VWQq 
Vq’ Vq 


Wq 


and we have 


relation which is also immediately verified for the case m= Mf 
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4. Bounds for the Determinants, depending on 2. and 4, 


\ determinant 


Mi ws mm, 


will be called an M/-determinant if all diagonal elements a, are positive, all elements off the 
main diagonal Mm,» (uv) are not positive and the determinant ./, as well as all principal 
coaxial) minors of .\/ of all degrees, are positive 

In what follows we will have to use a theorem given by the author in 1937 [1], and which 
will be formulated as the following: 

Lemma | If we have for the M-determinant (50) and a determinant H 


of order n, the inequalitie 8 


h = (ty, h., =m,, 


then 10. and we have 
H 
Tf (1), (M), re spectively, denote the matrices of the determinants H and M, the inverse matria 
of I] 18 majorized by the inverse matrix of (M) 


(H)"<«(M (32° 


Suppose now that in the determinant Q, given by (17) all a, are positive and m,, M, non- 
negative. Then it follows from (18) that Q, is a monotonically decreasing function of all m, 
and from (20) that M is also monotonically decreasing in all M,. Suppose now that for a 
certain set of values of a,, m,, M, the determinant 2,#0. Then replace the m, and 


V, corresponding to certain rows (»,, . . .,») by zeros; the determinant Q, cannot decrease 


and remains therefore positive; but then 2, becomes equal to the product of aa, . . . a, 
with the principal minor complementary to the set of indices »,, . . .,.».. Therefore, all principal 
minors of Q, are positive. We see that Q, is an M-determinant, if the conditions 


a,>0, m,20, M,20 (y a oe 


are satisfied and 2,#0. But now we can easily deduce the following theorem. 
A. Consider the set of all determinants A=\a,,| satisfying the conditions 


a,|sm, (v<y), a,|sM, (v>up), n), (54) 


where a, are n given positive constants and m,, M, are 2n—2 given nonnegative constants. Then 
in order that no determinant A of this set vanishes, it is necessary and sufficient thatQ,>0. If 


this condition is satisfied, we have 
A|2=Q,. (55) 


Proor: Since in the case 2, >0, Q, is an M-determinant, the sufficiency of our condition 


follows immediately from the inequality (52) mentioned above. 
If 2,—0, we can take obviously a,,—a, and a,,—=—m, or a,,——M, according as 
y<y or » >» and obtain a vanishing determinant A satisfying the condition (54). 
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Suppose now 2,<0, then it follows from the form (17) of Q, that @, becomes positive if 
all m, are replaced by zeros. There exists, therefore, such a positive f<1 that @, vanishes if 
all m, in (17) are replaced by tm,, but then we obtain a vanishing determinant A of our set 
in taking a,,—a, and a,, tm, or a@,, V/,, according as vy or y>u. The theorem A 
Is proved 

In specializing the matrix of Q, to A, and in assuming that in particular 0< m< MM, we 
obtain immediately from (22) and the theorem A 

B. Consider for two positive constants m, M, m< M, the set of all determinants A a,, of the 


nth order for which 
a - (u sn), - a,isM (: 56 


then, in order that no determinant A of this set ; ish. 4 s necessary and sufficient that for 
Ma VJ the neg uality 


Ml 


VW)’ a 


holds. and if this nequality is satisfied, we have for each determinant A of the set 


aici m)"(1+M)’ Ul " 
: VV UL ] V/ Me 


Tf m Vl. the condition hecomes 


and we ha ‘ 57 is vatished. 
Al=iA Lymi(1 m)’ 
From theorem B we can deduce the following theorem 


an) be a matria satisfying the conditions 


where m V/ are two constants with 0 m 


V/m Mm V/ 


oll, V/ VV ™ 


The n all fur damental roots of the matria A are contained in the set of the n closed circles described 
around the elements Gon with the radius é6(m.M The value (60) of 6(m.Mf) cannot be mproved 


if A=A 
Proor: Let X be a fundamental root of A so that the matrix \#—A 1s singular Put 


min A—4a,, a; we then have to prove that a<é(m,M If a=0, there is nothing to prove 


Suppose a >0, and consider the matrix 


For this matrix we have 











Therefore, we have by theorem B 


and therefore, 


M ™ m 


as m=" — \f=tin 5(m, M) 


5. Bounds of the Matrix A;' 


- 





For an n dimensional vector §=(x,, . . . .t,) the Hélder norm corresponding to the 


exponent p21 is given by 


4 
Ss 3 


We will only use the three cases corresponding to p=1, 2, « 


£ ryi- T \Zel, i. max J, , 


‘ 


We have among these three norms the following inequalities 


t] <itl <nlt 
so SS 8S i= Ss 

e . 
g S Ele Synig 
E, Synge Syng, 


Eee wvinlvrt t |_|" (p21) 


tr 
“ 
— 


(6) 


(63.8 ) 


(63.b) 


(63.¢) 


which are immediately verified. The left-hand inequality (63, ¢) implies the well-known 
inequality between the arithmetical mean and the arithmetical mean of the squares. If 


A= (a,,) is an nn matrix, we define its norm corresponding to the exponent p (lS pS @) by 


At 
A|,=max —" 
&0 Siz 
and denote it by .A|,. We will use here too only the cases p 


(64 


In applying to the definition (64) the formulas (63, b), and (63, c) we obtain immediately 


l 
113 Ales yn | 
yi 
l 
A = Ales niA 
yn 
For p=1, © the expressions of |A),, |A\, are easy to write down; 
very easy to prove, 
A|\,=max >> /4,, 
A|,.=max >} |a,, 


“ t 
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(65, a 


(65. b 


as is well known and 


(66, a) 


(66, b 











As to A», its expression is irrational; A) is the square root of the maximum fundamental 


root of the symmetric and nonnegative matrix AA*. Since the direct computation of |A|, is 


difficult in most cases, we prove the following estimate for A), 


Lema I] We have for any mati ix A 
max (|A/,,|A/.)S/AleSwy/Al,\Al., (67) 
” 
The first part of (67) follows from (65,a) and (65,b To prove the second part, we introduce 
the notations 
Su » > Qy|, b, + > ys g.e=lI, on 


The sum of the moduli of all elements in the uth row of AA* can be estimated as follows 


Dy Bua Bva| SD) Dy | Epc Ove| = > |All $8,|Al, S|Al.|A 


va 


Equation (67) follows now from the theorem of Frobenius that the modulus of each funda- 
mental root of a square matrix does not exceed the greatest sum of the moduli of the elements 


of this matrix in different rows 


If we now apply these results to the matrix A>' discussed in the section 3, we obtain from 


48), (49), (66,a) and (66,b) 


and from (67) and (65,a 
( a (.M +1)’ 
4 Satis . (69) 
yn(M+1)* M+1 A, 
In combining these inequalities with the results given in section 4, we obtain the following 


theorem 


1D Let A a, he a square matrix. of order n sctistying conditions (56 ‘ and let (57 he 
vsatished Then we hare 
(M+1 
| ~ A p=1,2,¢ 
and therefore for any vector = 
A, 
> t ) a 4 
Ag, (M+1)*-!'8" (p= 1,2,¢ (0 


Proor: It follows from (52°) of lemma I at once in virtue of (66,a) and (66,b) that 
, A 
: > 6 M+1 


(A~')* is majorized by A;'(A;')*. Py a well-known theorem of 


further, the matrix A 
A~')* is majorized 


Probenius, therefore, the maximum modulus of a fundamental root of A 


by that of A>'(A>')*, and so we also have 


By definition (64), therefore, in putting Ag=—n, we have 
4 

A 4 ! 

Vem 74!" 


and this is equivalent to (70 
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6. Bounus of the Matrix A, '— A.” 


We denote by A.” the matrix obtained from (21) in replacing m by zero: 








@ V/ M M) 
0 | M VV 
A (71) 
L0 0 0 1 J 





We will prove the following theorem 


EK. If 0<m Vl and (457) is satished, then the matrix A, a= is a nonnegative matriz 


and we have 


l m)" 4. = 
A A° max (1 +.M) paz ] i2 
Proor: Since A. satisfies the inequalities (51), it follows from (52°) that A, A*° is 
nonnegative. Denote by s{° and ¢{ the sums of the elements in the uth row and in the rth 
column of AS ~', and by s, and ¢, the corresponding expressions for the matrix A, A® 
We have 
:, i t,=t,—t p,y=l n 73 
It follows then from the formulas (48) and (49) that s, runs for w=n, n—1, ._ 1 through the 
same set of values as ¢, for yp=1, . . ., n. 
Formula (49) can be written in the form 
q V/ ] l m | VU 
1 74) 
4, 
Since A, becomes | for m=0, we have (72) and E is proved 
Discussion of A A We prove first that ¢ t, increases with vy as long as the 
condition 
mm A — 
> io 
VU 
is satisfied Indeed, put 
/ ] uM A. . 
j M l i6 
f ky | n A, 
In solving the three inequalities 
b,=1 (77 
with respect to A.|, we obtain correspondingly 
4, |= (1 77) ' 78 
. CG, - Mm \é 
> M 
The inequalities k,>1 (»—1, . n) are obviously satisfied in virtue of (75) as long as 
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. ven—l We have, therefore, in this case, in using (72) 
} Max (t,—1:%)=(1+ Myr) $4") _ gg ag : 
. A, l 6 
] A 6 / m 
: A A‘ 1+ M)-! ) : ( 
_ 1+ M i (p=1,0; |A,|<1 7) 79 


Condition (75) can be written by using the notation (27), as 6=m/.\/, and in virtue of (29) this is 
certainly satisfied if we have M@M,2 1, 


V+) nM—22>0 80 

Condition (80) becomes for n 3: W>0.5321 Forn>4 SU) sin any case satished. if we hare 
1.5 

V - n=a4 S| 
n 


Indeed, if the relation (80) is satisfied for a positive .\/, it is satisfied for any greater value 
because the coefficients of all positive powers of .V/ in the left-hand expression are greater than 


, or equal to 0. To prove the sufficiency of (81), it is sufficient to prove that 1.5/n)"23.5 
n=4 But here the left-hand expression is monotonically increasing with n, and this 
inequality follows, therefore, from (5.5/4)*=3.75 3.5 
Suppose now that we have 
m . , — 7 
l \/ A l 6 VW S2 


Then we have, since 0< m< 1/(n—1 


n l 
l m | 1.72 0nm 0< A< 1 
l m A 1.72 0nm l A, 
and therefore by (S82 l m -\|A,) <1.720nm+-m/M 
For n=>4 we have in the case (82), since the inequality (S81) is not satisfied 
] : n 
: VU L.d 
and, therefore, finally 
| ] "7 A 2? 39Onn ? ant 
On the other hand, since W<1.5/n, (1 \V/ <(] L5/ny"<¢ 1.481689, and, therefore 
from (72) and (73) 
/ | ) = hi | A hi ) 
/ ” <10.72n l : n>o p 
O< arn ( " 
" A AM 10.72 Mi ( , 4 Me :) Qf 
. (all 5 . n » 
{ "ie P V/ 


To obtain a lower bound, we take im (74) : ] We obtain, since for n >3, 0 A 








and, therefore, . 


— m ‘ 
A A »20.75n a p=1,°;n24). (85 
To obtain the exact value of A; AY ~', if (75) is not satisfied, we return to the inequali- 


ties (78) which are equivalent to (77 
Denote by no the smallest integer between 1 and n, such that we have 


0 > iA, CG, 1S Sn) (S6 | 
. 


The parts of this inequality implying oo or ¢, must be disregarded, that is to say, this inequality 


reduces to A o, for n 1 and to a, 2 (A, for no=n Then we see at once that 
max (t,—t? t..— te and therefore 
"A Ae P ‘ 1+ M)"e-'(1+m o—\A ' 
” ( ) » ) R7 
o mq / 
For m= M we have no=1, and therefore by (23) 
n—l1)m ae 
._ m=M< (87° ¥ 
4, , l (ni l)m ( . f n 1 ) ) 


7. Bounds of A~'— A” 


’ 
In order to obtain the theorem corresponding to [f for a determinant A satisfying condi- 
tions (56), (57) of theorem B, we prove first the following important lemma, which generalizes 
considerably the relation (52°) of lemma I 
Lemma III (Consider n positive numbers a), , a, and 2n? 2n nonnegative numbers 
Myr, Exe (FV; pwv=—l, n), such that the matrices 
f a m € m € % 4 
m € Qe m € 
i SS 








Vl (89 








a Mm nt M »2 a, : i 


are \f-matrices. | 
(Consider n constants A, vi = i) a , such that 


A,| =a (n= 1, ~n 90 


| 
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a 


and 2n?—2n constants a,,», 6,, (ux¥tv: ww. _n) such that 


ys =™M,,», b, Gun! Sen (ux: uy l, nh). (91) 


and form the trvo matrices 





[A h b,, ) 
b Ay b 
| (9°? 
Vb b.» Aus, 
(A a a a 
a Ay a 
A* (O:3) 
Ldn, Qn A J 





Then we ha ¢ 
A A-\e M Wwe (94) 
Proor: We write M—=P—T, M* P—T”, where the matrices 7, 7 have zeros 


in the main diagonal and off the main diagonal, respectively, the nonnegative elements m,,-+-« 


m,, and where P is the diagonal matrix 


fa \ 








\ AanZ 


We can develop the inverse of P— T=P(£—P-'T) (F is the unit matrix) in the following was 
P—T zur aye 96 


The convergence of this development and the validity of (96) follow easily from the fact 
that the determinant of the matrix P—tT does not. vanish for t_<1, which follows immediately 


from lemma I. The corresponding development holds also for (/?~ 7 and we obtain 
therefore - 
M Vl DP OT) —(P 7 / 
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The elements of P T are here (Why €yy) ay OT ZeTOs, and those of r T" are W,, @, Or Zeros: . 


therefore all elements of the matrices 
r7)’ pT; a=1,2. Q8 


are polynomials in m,,/a, and «¢,,/a, with nonnegative coefficients. Denote by ( the diagonal 


matrix 








1 () . 
J) 
‘ 
) 
V A.J 
nd write A=V—S A —S, where the elements of S are 4,, or zeros, and those of 
SS , a or zeros snce therefore we have 


we have 
. . " 
j S*(0-'S)\*Q A > -(0)-*4 ij 
| A S*1(0-'S) O-'S8 0 90 
—_ ‘ ‘ ‘ 
But now the elements of (Q-'S)* -'S®)* are obtained from those of (98) in substituting 
there a,,/A, instead of m,,/a, and (b,,—a,,)/A, instead of €,,/a Further, we have by (90 
and (91 
a m, b,, a € 
a a, 1, a 
and therefore, since the coefficients in (98) are as already mentioned not negative, 
v's g-'S©) a 8 Pz 
and | 
S*1(Q0-'S)* o's VeS‘{(P-'T) P-'7 “ 100) 


But from (100) by virtue of (90) it follows that the development (99) is majorized by (97) 
and our lemma is proved 
Under the conditions of theorem B in section 4, if the inequalities (56) and (57) are satisfiec, 
we can applv lemma III in replacing the matrix .V by A,, 17° bv A® We obtain 
j EF. Lnder the conditions of theorem B. if the inequalities 56) and (57) are satished, we hare 
i l m)"~* A, 
| A A s\A A p= max(1+ MVM p= 1,< 101 
nites A, 
where the values and the estimates for the right-hand EXPressit n in (101) are obtained from the 
formulas (79) te (87° 


Lemma III can be applied to many problems similar to that solved by theorem Ff. For 


+ rereer ener 
SEO 


instance, in the theory of the solutions of linear equations, the following problem has to be 
dealt with, although its complete discussion is usually avoided q 
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es 
ae 


Prat aed 





Consider a ‘triangular’ system of linear equations 


> Osho Pp p=! " 


where the coefficients yy al 


‘only approximate values to the ‘“‘true”’ 


we have generally 


iow far is the solution of the system 


influenced if the 6,, are replaced by a,,' 


102 


values 5, Suppose that 


108 


104 


If we denote the matrix of (102) by A and that of (104) by A, the question can be 


answered by giving estimates of A A“ p=l,¢ Suppose that we have generally 


a l. a,,| s M, u>v; ob | ! 


and consider the matrix 





) ) ] 


106 





J: 


where all elements to the left of the main diagonal are zeros, and all elements to the right of 


the main diagonal in the wth row are equal to VW, 
We obtain from lemma II] at once the majorization 


A A <A(M,+«. , M,_st+e A(M,, 
To obtain the inverse of the matrix ACV, . M,_,), consider the 
J V (a J u a= | 


. > (mu l n 
1+ M,=N,, N,N 4 N =F, (p=1, , n—!) AN 
Then (107) becomes 
r.=M,3s,.,+4 (u==1 n 
= N Sei t+ Uy, w=) n 
where s 0. Dividing (111) by P,, we obtain 
x x ' u 


” m1 


P +1 
fn Set Fo 
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system of linear equations 


n 107) 


(10S 
l, | 109) 








and therefore 














“—" 
. >>. Me. 
a Podu P. 
introducing this in (110) we have 
r,=u,+M,P, > p w=, 2; P4,;=0) 112) 
We obtain therefore for the inverse of our ACW, _ M,_,), if & denotes the unit matrix 
ACM,, , M,.1) E+DTrT. 113 
where J is the diagonal matrix 
(M,P, 7 
VP 
0) 
dD ; 114 
0) 
MF 
\ OJ, 
and 7 the triangular matrix 
i i= ca 
ro F 
Q l ] 
P P 
0 , 
ij P 115 
0) 
3 0 7 
The expressions 
A(M, +e, , M,_site A( VM — 116 
obtained from (113), (114), (115) are rather unwieldy; however, we shall obtain for its norm 


corresponding to p=1 (ef. 66,a) a very simple and elegant expression. 
Indeed, if we take the sum of the elements in the rth column of ACM, . , M,_,)“' and 
denote it by t’, we have 


' ee _— 
14 p S MPa 1+ > MM, (dl + M,). 
If we now write out ¢.,, we obtain 


’ 1 ’ 
Lo i+{ SM, i 


p=l c=p+l 


(1-4 M. | (1+M,)+M,; 
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ete tenis sak ade 











ere womens eres 
—_s—- 


and in comparing this with the expression for ¢t! we see that we have t,, 1+ \/,)t! and 
therefore 


(= N(1+M, 117) 


We obtain now for the sum of the elements in the rth column of (116 


In multiplying this by 1+.M, ve obtain 


“a 
< 


and comparing this with ¢, 


t.s(1+3+-06+¢ 004M, 118) 


Therefore, t, is monotonically increasing with v, and we see that the norm of (116) cor- 
responding to p=1 has the value 
i (1+ M,+e (1 +M, 119 
and obtain therefore 


A“—A < I (1 |M,+e)— N(1+M 120 


as the solution of our problem. In applications it may be better to use the recurrence formula 
118) If all \/, have the same value VV, the expression (113) coincides with that obtained in 
section 5 for A-' for m=0. But in this case we see from (48) and (49) that the row sums run 
through the same values as the sums in the columns. We obtain, therefore, in this case the 


expression 
V/+-1+-e)” VJ+-1)’ hi Le( MM ] He). i) “ is 


as the norm of (116), both for p=) and p 


8. Linear Systems With a Nearly Triangular Matrix 


The results of the preceding sections give the means to discuss the following problem 
I | a) 


concerning the system (1) in the introduction under the conditions (2) and the “triangular” 
svstem 
Sayed, UP p=] , >, 12] 

with the matrix A“ In discussing this problem we can obviously assume that 
é,,.=1 (sa=!1, n Then the difference between the solution of (1) and that of (121) is 
viven by the vector (A A"), 9 Y, Yn), and the norm of this vector corresponding 
to one of the indices p=1, © does not exceed A A n\,» (p=1,@), and can indeed for 
suitable choice of the vector » attain this limit. Therefore, the norms |A A measure 


the error committed in replacing the system (1) by (121 
For an e>0.M~>0 being given. how small must m™>0O be taken in order that we have 


A A‘ Se p=1,)? (122) 








2) (2: 


> 
ae), 


If we introduce the quantities |A, and 6 corresponding by (2 


\/, we obtain from (79) and (94) the condition 


. € 5 UL 
Saupe, (ll 51-97) 


as long as the condition (75) is satisfied and therefore certainly as long as 


On the other hand, we have by (30) 


6 


M.+6 M 


and from (123) and (124 


€ 


M(1+ M)""'+(M, +57). 


We see that it will be sufficient to take 


msm ty (mls) 


to solve our problem, for instance, if we have M21.5/n (n24) or M20.5321 


and 


V/ 


(97 


~f 


N 


to m and 


n(n2az4 


small values of « obviously, only the first term in the denominator is essential, and we have 


i "(nA “(n 
Kin, Me, Kin,M VATS TAT 


re 


~/ 


Tables 1 and 2 give the values of K(n,M) for a set of integer Mf from 1 to 10 and some values 


of V/ ] 5 hi We have obv iously 
l 


K(n,M) <5 


TARLE 1 Kin, M 


V n ; , j n ) n 6 n | n bed n 
! 0625 Oll4 2)240 2)548 3)130 1)316 5 
2 vlll 2)103 3) 106 $)115 >) 126 6)140 7 
3 2)347 3) 193 4)116 6)719 7)448 8) 280 9 
{ 2)143 $)526 5) 206 7)821 8)328 9)131 11 
y 604 1) 182 6)498 7)138 9) 383 10) 10¢€ 12 
“i 3)378 5)736 6)149 8)?04 19)619 11)126 13 
7 3)223 5)336 7)522 O)R815 10) 127 12) 199 14 
s 3)140 5) 168 7 )207 9)255 11)315 13 )389 15 
a 1)926 6)903 8 )900 19)900 12)900 14)900 16 

10 1)636 H)515 &)424 10)351 12)290 14)239 16 
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}900 
198 


(128 


10 


5)192 
8)172 
10) 109 
12)210 
14)821 


15)526 
16) 486 
17) 592 
18)900 
18) 164 





a raRLe 2 K (in, M 


V ; , j , ) n=6 , 7 x 9 10 
l 
“ 2 032s O107 0121 

7 O66 0399 01906 0101 2)533 2)28%) 

: 117 0438 OLSO 9)788 2)359 2)169 800 

) 254 0718 0241 2)801 2)349 2)142 3) 502 3)252 

6 18] 0464 0141 2)470 2)165 3) 602 $)224 1)848 

he 7 134 0313 OO864 2) 260 3)825 $)270 1)002 1)205 
S 102 0218 00548 2)150 3)431 $)127 1)382 1)116 

q 0789 0156 00359 3)804 3)233 1)623 1) 169 >) 462 

0 vV625 O1l4 00240 3)548 25120 1316 5)778 102 

15 n=—20 n=25 n — 30 n= 35 , 0 , 5 ’ 0 

l O157 2)439 2)138 3)469 3) 166 1) 604 $)224 »)S41 

2 2)137 3) 188 $)281 5) 439 6)697 6)112 7) 180 8) 200 

; 3) 167 $)112 6)793 7)570 8)413 9) 209 10)217 LL)157 
} $)242 6)809 7)277 9)957 10)331 11)114 13)305 14)137 
9)399 7)681 S)118 10) 204 12)354 14)613 15) 106 17) 184 
( 6)729 8 )658 10) 597 12)543 14)494 16)449 18) 400 20)372 
7 6)146 9)721 11)357 13)177 16)879 18)436 20) 216 22)107 
8 7)317 10) 886 12)248 15)605 17) 195 20)545 22)153 25)428 
9 £)743 10) 121 13) 197 16)322 19)525 22) 856 24)140 27 ) 228 
4 1.0 8) Se 11)182 14)178 17)173 20) 169 23)165 26) 162 20) 158 


The bound in (126) is obviously the “best’’ under the condition (75), save that the factor 
1/\/ of « in the denominator could be replaced by an (unknown) fraction of it 
If \f<1.5/n, the use of the general formula (87) is very cumbersome. We can, however 


obtain a good working limit for m in the following way. If the inequalities (2) are valid for 
an .WM<1.5/n=—M™, then they are also satisfied if \/ is replaced by MM, but then the limit 
f 


for m obtained for \/” is also sufficient for our V/ We obtain therefore in this case the 
sufficient condition for (122) in the form 


} saith l 
| Mo (1+M +(Mir+yy)e 
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Table 3 vives the Vaiues of K n. V/ . for n L.d- . 50 


TARLE 3 V Wr L.5/n 
n l VJ nM—1 Kin, M n I W)"*—nM—1 Kin, M 
3 S71 255 27 1. 805 00755 
i 1. O74 134 8 1. Sil 00723 
5 1. 213 O866 at) 1. 817 00694 
6 1.315 0623 30 1. 822 00667 
7 1. 393 0480 31 1. 827 00642 
Ss 1. 454 O387 32 1. S31 OO619 
9 1. 504 0323 33 1. 836 00597 
10 1. 546 0276 34 1. 840 00577 
11 1. 580 0240 35 1. 844 00558 
12 1. 610 0213 36 1. 848 00540 
13 1. 635 0190 37 1. 851 00524 
14 1. 658 O172 38 1. 854 00508 
15 1. 677 O1L5ST 39 1. 857 00493 
16 1. 695 0144 10 1. 860 00480 
17 1. 710 0133 +1 1. 863 00466 
18 1. 724 0124 42 1. 866 00454 
19 1. 736 O116 13 1. 869 00442 
20 1. 748 0109 44 1. 871 00431 
21 1. 758 0102 5 1. 873 00420 
22 1. 768 00965 16 1. 876 00410 
23 1. 776 00914 17 1. 878 00401 
24 1. 784 OO869 18 1. 880 00391 
25 1. 792 00827 49 1. 882 00383 
26 1. 799 00789 50 1. 884 00374 


The expression for m, can be written in introducing the value of M\ as 


(14h ((14 48) —2.5)4((1419)—-1.5)e 


We observe now that for any positive a and positive z the expression (1+ 2/n)""* mono- 
tonically increases and tends to ¢ if the positive n increases monotonically to ~. Indeed, if 
we put u=1/n, take the logarithm of this expression, differentiate it with respect to u, and 
multiply by u?, we obtain 

r(u—au*) 


—log (1 ur) 
l+ur 


but this expression vanishes for u=0 and decreases for positive u, since its derivative ts 


ru 


We see, therefore, that 


and our bound for m can therefore be replaced by 


» 


» € € 


2n e'*(e' 2.5)+(e'°—1.5)e~ (5.934+1.99€)n 
we can replace therefore the condition m <m, by the simpler condition 


ms™Mb, : . . (130) 
7 


1 6+-2e 
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from (87 


This is a “‘best’’ condition, whet 
however, in comparing (130) with 
factor greater than 6 


From the condition (130) we 


Indeed, if POSILIVe Mo: 1/2n 1s given, 


On the other hand, we obtain at once the solution of our problem in the case m 





l € l “ 
; ( m M ) 131 


eas the condition (130) can still be improved. We see 
131), that the bound in (130) cannot be improved by a 


‘an finally derive the inequality (7) of the introduction 


we can solve (130) with respect to e, 


onm 


l 2n Mm 


and obtain therefore, in applying (130) for m=mp, the inequality 


A 


which holds for all positive mp<1/2n 


fl] A. Ostrowski, Uber die Determinant 


69—96 (1937 


onm, 
A* 


1—2nmy 
and gives the inequality (7) if we replace here mo by m 
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